Interaction-induced magnetoresistance in a two-dimensional electron gas 
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We study the interaction-induced quantum correction 5a a p to the conductivity tensor of electrons 
in two dimensions for arbitrary Tt (where T is the temperature and r the transport scattering time), 
magnetic field, and type of disorder. A general theory is developed, allowing us to express 8a a /3 
in terms of classical propagators ("ballistic diffusons"). The formalism is used to calculate the 
interaction contribution to the longitudinal and the Hall resistivities in a transverse magnetic field 
in the whole range of temperature from the diffusive (TV <C 1) to the ballistic (TV > 1) regime, both 
in smooth disorder and in the presence of short-range scatterers. Further, we apply the formalism 
to anisotropic systems and demonstrate that the interaction induces novel quantum oscillations in 
the resistivity of lateral superlattices. 

PACS numbers: 72.10.-d, 73.23.Ad, 71.10.-w, 73.43.Qt 



I. INTRODUCTION 

The magnetoresistance (MR) in a transverse field B 
is one of the most frequently studied characteristics of 
the two-dimensional (2D) electron gasi 2 -. Within the 
Drude-Boltzmann theory, the longitudinal resistivity of 
an isotropic degenerate system is .B-independcnt, 



Pxx(B) = p = (c 



(1.1) 



where v is the density of states per spin direction, vf the 
Fermi velocity, and r the transport scattering time. De- 
viations from the constant p X x{B) are customarily called 
a positive or negative MR, depending on the sign of the 
deviation. There are several distinct sources of a non- 
trivial MR, which reflect the rich physics of the magne- 
totransport in 2D systems. 

First of all, it has been recognized recently that even 
within the quasiclassical theory memory effects may lead 
to strong MR 3 i 4 i 5 i 6 i 7 i 8 i 9 . The essence of such effects is 
that a particle "keeps memory" about the presence (or 
absence) of a scatterer in a spatial region which it has 
already visited. As a result, if the particle returns back, 
the new scattering event is correlated with the original 
one, yielding a correction to the resistivity i|l.l|) . Since 
the magnetic field enhances the return probability, the 
correction turns out to be B-dependent. As a prominent 
example, memory effects in magnetotransport of com- 
posite fermions subject to an effective smooth random 
magnetic field explain a positive MR around half-filling 
of the lowest Landau level 7 -. Another type of memory 
effects taking place in systems with rare strong scatter- 
ers is responsible for a negative MR in disordered antidot 
arrays 3 * 4 *^^. However, such effects turn out to be of a 
relatively minor importance for the low-field quasiclassi- 
cal magnetotransport in semiconductor heterostructures 
with typical experimental parameters, while at higher B 
they are obscured by the development of the Shubnikov- 
de Haas oscillation (SdHO). 

Second, the negative MR induced by the suppression 
of the quantum interference by the magnetic field is a 



famous manifestation of weak localization 1 . While the 
weak-localization correction to conductivity is also re- 
lated to the return probability, it has (contrary to the 
quasiclassical memory effects) an intrinsically quantum 
character, since it is governed by quantum interference of 
time-reversed paths. As a result it is suppressed already 
by a classically negligible magnetic field, which changes 
relative phases of the two paths. Consequently, the cor- 
responding correction to p xx in high-mobility structures 
is very small and restricted to the range of very weak 
magnetic fields. 

Finally, another quantum correction to MR is induced 
by the electron-electron interaction. While this effect is 
similar to those discussed above in its connection with the 
return probability (see Sec. II VI below), it is distinctly dif- 
ferent in several crucial aspects. In contrast to the mem- 
ory effects, this contribution is of quantum nature and is 
therefore strongly T-dependent at low temperatures. On 
the other hand, contrary to the weak localization, the in- 
teraction correction to conductivity is not destroyed by a 
strong magnetic field. As a result, it induces an apprecia- 
ble MR in the range of classically strong magnetic fields. 
This effect will be the subject of the present paper. 

It was discovered by Altshuler and Arono\i that the 
Coulomb interaction enhanced by the diffusive motion of 
electrons gives rise to a quantum correction to conduc- 
tivity, which has in 2D the form (we set ks = h = 1) 



e 

2^2 



1- -JF ^Tt, 



Tt < 1. 



(1.2) 



The first term in the factor (1 — |jF) originates from 
the exchange contribution, and the second one from the 
Hartree contribution. In the weak-interaction regime, 
k <C hp, where k = Aite 2 v is the inverse screen- 
ing length, the Hartree contribution is small, T ~ 
(k/^f) ln(fci?/K) <C 1. The conductivity correction ljl.2|l 
is then dominated by the exchange term and is negative. 
The condition Tt <C 1 under which Eq. I|1.2fl is derived^ 
implies that electrons move diffusively on the time scale 
1 /T and is termed the "diffusive regime" . Subsequent 



2 



works io i ii showed that Eq. |T2| remains valid in a strong 
magnetic field, leading (in combination with Sa xy = 0) 
to a parabolic interaction-induced quantum MR, 

SPxJB) ( 3 \ (uo c t) 2 - 1 , 
F 1 - ~ ( 1 — -T I , — lnTr, Tr < 1, 
Po V 2 / 

(1.3) 

where w c = eB/mc is the cyclotron frequency and 
I = vpT the transport mean free path. In- 
deed, a T-dependent negative MR was observed in 
experimentsi2ii&±!ii5i±& and attributed to the interac- 
tion effect. However, the majority of experimentsi2iiii4 
cannot be directly compared with the theoryiiS*ii since 
they were performed at higher temperatures, Tr > 1. 
(In high-mobility GaAs heterostructures conventionally 
used in MR experiments, 1/r is typically ~ 100 mK and 
becomes even smaller with improving quality of sam- 
ples.) In order to explain the experimentally observed 
T-dependent ne gative MR in this temperature range the 
authors of Refs. Il2fl3l conjectured various ad hoc exten- 
sions of Eq. (|1.3(l to higher T. Specifically, Ref. 12 conjec- 
tures that the logarithmic behavior (jl.3|l with r replaced 
by the quantum time r s is valid up to T ~ 1/t s , while 
Ref. proposes to replace In Tr by — 7r 2 /2Tr. These 
proposals, however, were not supported by theoretical 
calculations. There is thus a clear need for a theory of 
the MR in the ballistic regime, T > 1/r. 

In fact, the effect of interaction on the conductiv- 
ity at T > 1/t has been already considered in the 
literaturo^i 18 ! 19 ! 2 "! 21 ! 22 ! 23 . Gold and DolgopoloOi an- 
alyzed the correction to conductivity arising from the 
T-dependent screening of the impurity potential. They 
obtained a linear-in-T correction 5a ~ e 2 Tr. In the 
last few years, this effect attracted a great deal of in- 
terest in a context of low-density 2D systems showing a 
seemingly metallic behavios^Sii, dp/dT > 0. Recently, 
Zala, Narozhny, and AleinesiS^i^i^ developed a system- 
atic theory of the interaction corrections valid for arbi- 
trary Tr. They showed that the temperature-dependent 
screening of Ref. ^3 nas m f ac t a common physical origin 
with the Altshuler-Aronov effect but that the calculation 
of Ref. took only the Hartree term into account and 
missed the exchange contribution. In the balli stic range 
of temperatures, the theory of Refs. Il8ll9l20l predicts, 
in addition to the linear-in-T correction to conductivity 
o xx , a 1/T correction to the Hall coefficienljiS p X y/B at 
B — ► 0, and describes the MR in a parallel fieldSS. 

The consideration of Ref. is restricted, how- 

ever, to classically weak transverse fields, lo c t <C 1, and 
to the white-noise disorder. The latter assumption is be- 
lieved to be justified for Si-based and some (those with 
a ver y large sp acer) GaAs structures, and the results of 
Refs. liMlflll2(l have been by and large confirmed by most 
recent experiments^SLSSiSiLSiiSiii^ on such systems. On 
the other hand, the random potential in typical GaAs 
heterostructures is due to remote donors and has a long- 
range character. Thus, the impurity scattering is pre- 
dominantly of a small-angle nature and is characterized 
by two relaxation times, the transport time r and the 



single-particle (quantum) time r s governing damping of 
SdHO, with r ^> r s . Therefore, a description of the MR 
in such systems requires a more general theory valid also 
in the range of strong magnetic fields and for smooth 
disorder. [A related problem of the tunneling density of 
states in this situation was studied in Ref. l33l ] 

In this paper, we develop a general theory of the 
interaction-induced corrections to the conductivity ten- 
sor of 2D electrons valid for arbitrary temperatures, 
transverse magnetic fields, and range of random poten- 
tial. We further apply it to the problem of magneto- 
transport in a smooth disorder at lq c t 1. In the bal- 
listic limit, Tr 3> 1 (where the character of disorder is 
crucially important), we show that while the correction 
to p xx is exponentially suppressed for oj c -C T, a MR 
arises at stronger B where it scales as B 2 T~ 1 / 2 . We also 
study the temperature-dependent correction to the Hall 
resistivity and show that it scales as BT 1 / 2 in the bal- 
listic regime and for strong B. We further investigate a 
"mixed-disorder" model, with both short-range and long- 
range impurities present. We find that a sufficient con- 
centration of short-range scatterers strongly enhances the 
MR in the ballistic regime. 

The outline of the paper is as follows. In Sec. [H] 
we present our formalism and derive a general formula 
for the conductivity correction. We further demonstrate 
(Sec. Ill C|l that in the corresponding limiting cases our 
theory reproduces all previously known results for the 
interaction correction. In Sec. II I II we apply our formal- 
ism to the problem of interaction-induced MR in strong 
magnetic fields and smooth disorder. Section IIVI is de- 
voted to a physical interpretation of our results in terms 
of a classical return probability. In Sections Ivl and IVll we 
present several further applications of our theory. Specif- 
ically, we analyze the interaction effects in systems with 
short-range scatterers and in magnetotransport in mod- 
ulated systems (lateral super lattices). A summary of our 
results, a comparison with experiment, and a discussion 
of possible further developments is presented in Sec. lVIll 
Some of the results of the paper have been published in 
a brief form in the Letted. 



II. GENERAL FORMALISM 

A. Smooth disorder 

We consider a 2D electron gas (charge — e, mass m, 
density n e ) subject to a transverse magnetic field B and 
to a random potential u(r) characterized by a correlation 
function 

(u(r)u(r')) = tuflr - r'|) (2.1) 

with a spatial range d. The total (rT 1 ) and the transport 
(r _1 ) scattering rates induced by the random potential 
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are given by 



i = 2™ r *t W w, 

t s Jq 2tt 



= 2-kv 



2tt 



W(0)(1 -cos0), 



(2.2) 
(2.3) 



where W{4>) = w(2kF sin t) is the scattering cross- 
section. We begin by considering the case of smooth 
disorder, fcpd 3> 1, when r/r s ~ (kpd) 2 3> 1; generaliza- 
tion onto systems with arbitrary r/r s will be presented 
in Sec. Ill Bl We assume that the magnetic field is not 
too strong, lo c t s < 1, so that the Landau quantization is 
destroyed by disorder. Note that this assumption is not 
in conflict with a condition of classically strong magnetic 
fields (co c t 3> 1), which is a range of our main interest in 
the present paper. 

We consider two types of the electron-electron interac- 
tion potential Uo(r): (i) point-like interaction, Uo(r) = 
Vo, and (ii) Coulomb interaction, U~o(r) — e 2 /r. In or- 
der to find the interaction-induced correction 8a a p to the 
conductivity tensor, we make use of the "ballistic" gener- 
alization of the diffuson diagram technique of Ref. Q. We 
consider the exchange contribution first and will discuss 
the Hartree term later on. Within the Matsubara formal- 
ism, the conductivity is expressed via the Kubo formula 
through the current-current correlation function, 



milk 

1 r l l T 



-j£ j o dr I d'r(T Tja (r,T) ]fj (0,0))e* 



(2.4) 



where Qk = 2nkT is the bosonic Matsubara frequency. 
Diagrams for the leading-order interaction correction are 
shown in Fig.^and can be generated in the following way. 
First, there are two essentially different ways to insert an 
interaction line into the bubble formed by two electronic 
Green's function. Second, one puts signs of electronic 
Matsubara frequencies in all possible ways. On the third 
step, one connects lines with opposite signs of frequencies 
e„ > 0, e m < by impurity-line ladders (which are not 
allowed to cross each other). Finally, in the case of the 
diagram a, where four electronic lines form a "box" , one 
should include two additional diagrams, b and c, with an 
extra impurity line ("Hikami box"). 

The impurity-line ladders are denoted by shaded 
blocks in Fig. we term them "ballistic diffusons" . For- 
mally the ballistic diffuson is defined as an impurity aver- 
age (denoted below as (. . .)i mp ) of a product of a retarded 
and advanced Green's functions, 



X>(ie m ,ze n ;r 1 ,r 2 ;r 3 ,r 4 ) 
= 0(-e m e„)(G(ri, r 2 ; ie m )G(r 3 , r 4 ; ie n )) imp . 




FIG. 1: Exchange diagrams for the interaction correction to 
<j a p. The wavy (dashed) lines denote the interaction (impu- 
rity scattering), the shaded blocks are impurity ladders, and 
the +/— symbols denote the signs of the Matsubara frequen- 
cies. The diagrams obtained by a flip and/or by an exchange 
+ <-> — should also be included. "Inelastic" part of the dia- 
grams /, g is canceled by a contribution of the Coulomb-drag 
type, see Appendix 1X1 



Following the standard route of the quasiclassical formal- 
ism^iSSiS, we perform the Wigner transformation, 

T>(ie m - ie„;Ri,pi;R 2 ,P2) = / drdr' 

Xe -i[Pi-(e/c)A(Ri)]r e -i[p 3 -(e/c)A(R 2 )]r' 

xX>(ie m ,ie„;Ri,r;R 2 ,r'), (2.6) 

where Ri = (r 4 + ri)/2, R 2 = (r 2 + r 3 )/2, r = r 4 - ri, 
and r' = r 2 — r 3 . Note that the factors depending on the 
vector potential make the ballistic diffuson l|2.6|l gauge- 
invariant. Finally, we integrate out the absolute values 
of momenta pi j2 and get the final form of the ballistic 
diffuson 



P(«w;;Ri,ni;R 2 ,n 2 ) 
_ _J_ f Pidpi f p 2 dp 2 
~ 2ttv 



2tt 



2tt 



X>(w^;Ri,pi;R 2 ,p 2 ), 



(2.7) 



which describes the quasiclassical propagation of an elec- 
tron in the phase space from the point R 2 , n 2 to Ri, ni. 
Here n is the unit vector characterizing the direction of 
velocity on the Fermi surface. The ballistic diffuson sat- 
isfies the quasiclassical Liouvillc-Boltzmann equation 



9 

0q) + Uc-^r- 



c 



V(iuji,q] 4>, (j)') 



271-5(0-0'), (2. 



where (0 g ) is the polar angle of n (q) and C is the col- 
(2.5) lision integral, determined by the scattering cross-section 
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FIG. 2: Diagrams obtained by a flip and/or by an exchange 
+ <-> — from the diagram d. 



W(n, n'). For the case of a smooth disorder, the collision 
integral is given by 



C 



i d 2 



(2.9) 



In contrast to the diffusive regime, where V has a uni- 
versal and simple structure X>(zwz,q) = l/(Dq 2 + \coi\) 
determined by the diffusion constant D only, its form in 
the ballistic regime is much more complicated. We are 
able, however, to get a general expression for Sa a /3 in 
terms of the ballistic propagator T>(u>i, q; n, n'). 

The temperature range of main interest in the present 
paper is restricted by Tt s <C 1, since at higher T the 
MR will be small in the whole range of the quasiclassical 
transport lo c t s <C 1 (see below). In this case the ladders 
are dominated by contributions with many (3> 1) impu- 
rity lines. We will assume this situation when evaluating 
diagrams in the present subsection. A general case of 
arbitrary Tt s and r s /r will be addressed in Sec. IIIBl 

We start with the diagrams d and e that give rise to the 
logarithmic correction in the diffusive regimei. Let us fix 
the sign of the external frequency, fl k > 0. Each of the 
diagrams d and e generates four diagrams by a flip with 
respect to the horizontal line or by exchange + <->—, see 
Fig- El Consider first the diagram d ++ . There are two tri- 
angular boxes containing each a current vertex and three 
electron Green's functions (Fig. [3]) . In the quasiclassical 
regime lo c t s <C 1 one may neglect the effect of magnetic 
field on the Green's functions (keeping oj c in the ballis- 
tic propagators only). Furthermore, using Tt s <C 1, we 
neglect the difference in momenta and frequencies in the 
Green's functions, since typical values of frequencies ifl k , 
iu)i and momenta q carried by the ballistic diffusons are 
set by the temperature. Each triangle then reads 



r Q (n) = - 



pdp 



P n c 



m J 2tt + i/2T s ) 2 (-Z P - i/2r s ) 
i2ni / T 2 evFn a , (2-10) 



where £ p = p 2 /2m — [i. Combining the triangles with 
the three ballistic propagators separated by the impurity 
lines (see Fig.OJ, we obtain the following expression for 
the electronic part of the diagram 



(2™) 3 J H ^23(iw,, q; cf>, h)W(4>x - 2 )r Q ( 

i — 1 

xW((j} 2 - <j>s)T>{iuJi - ifl k , q; fo, 4>4)W(<f>4, - 5 
xi>(0 5 )w^ s ~ faWwi, q; ^6, 0') 



2) 



47TC7o 



Bipfau -^fc,q;<£, 4>')- 



(2.11) 



In what follows we will use for brevity a short-hand no- 
tation 

{2-xvfVWY a WVWT P WV 

for the l.h.s. of l|2.11|l and analogous notations for 
other structures of this type. Making use of the small- 
angle nature of scattering in a smooth random poten- 
tial, we can replace the W(4>i — <fij) factors in (|2.11|l by 

{vr^S^i - (j>j), yielding 

~ -T>{iu h q)n a r>(iw; + itt k , q)npT>(iu}i, q). 

In the exchange term (calculated in the present subsec- 
tion) this structure is further integrated over the angles 
<p and cj)' , 

fl^(iw,,in fc ,q) = (B*,(iw,,ifi fc ,q;(^)>. (2-12) 

The angular brackets (. . .) denote averaging over velocity 
directions, e.g. 



(n, r 'Dn„ 



#1 d(f> 2 
~2tt^2it 



cos <j>iD{w u q; 0i,0 2 



The fermionic frequencies obey the inequalities e m > 0, 
e m — lui < 0, and e m — £l k > 0, which implies uji > 
e m > f2fc, so that the summation over e m gives the factor 
(w, - n k )/2nT. 

The diagram d has the same structure (both tri- 
angles have opposite signs, thus the total sign remains 
unchanged), but the frequency summation is restricted 
by e m < 0, e m - wj > 0, and e m — Sl fc < 0, yielding the 
factor —ll>i/2itT in the conductivity correction. The di- 
agrams d' ++ and d! obtained from d ++ and d by a 

flip (or, equivalently, by reversing all arrows) double the 
result. Combining the four contributions and changing 

sign of the summation variable, cjj — > —lui in d and 

d'__ terms, we have 



87rcr T 2 f d 2 q 



U>1 



T 

n k 



2irT 



fl k J (2tt) 2 
U(iwi,q)B*p(-iwi,iSlk, q) 



(2.13) 
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FIG. 3: Diagram d drawn in a different way in order to visu- 
alize the structure of Eq. 12. lit . The dashed frame encloses 



the electronic part SJL 



where U{iu>, q) is the interaction potential equal to a con- 
stant Vq for point-like interaction and to 



U(ioJi,q) = 



1 



2v g + K[l-|wi|(23(q,iwi)>] 



(2.14) 



for screened Coulomb interaction. In l|2.13|) we used 
the fact that U(q, — iuii) — U(q,iuii) and £>(q, — iuJi) = 
P(q, iu>i). Equation H2.14(l is a statement of the random- 
phase approximation (RPA), with the polarization oper- 
ator given by 



n(*w I> q)=2i/[l-|wi|(I>(<wi > q))]. 



(2.15) 



The first term (unity) in square brackets in Ij2.15|) comes 

from the H — h and contributions to the polarization 

bubble, while the second term is generated by the H — 
contribution (ballistic diffuson). 

The diagrams e are evaluated in a similar way. In all 
four diagrams of this type one of the electron triangles 
is the same as in diagrams d while another one has an 
opposite sign. The structures arising after integrating out 
fast momenta in electron bubbles coincide with those of 
d-type (B^g). Summation over the fermionic frequency 
e m is constrained by the condition u>t > e m > ft k for all 
the e-type diagrams. The correction due to the diagrams 
e therefore reads 

x e rn \ 871-ctoT 2 f d 2 q ^ ui - 0, k 

^ ( ^ ) = ^^i(2^ D-w"^'^ 

x [B d ap {-iio h ifl k , q) + B^(iui,%a k , q)] . (2.16) 

We see that the first term in square brackets in (|2.16(l 
cancels the first term in <|2.13|l . Thus, the sum of the 
contributions of diagrams d and e takes the form 



T flk 



J2 WiS^Ctwj.tfifc) 

Wi=0 



(2.17) 



where we introduced a notation 
d 2 q 



(2tt) 2 



U(w h q)B£JiuJi,iQ, k ,q), 



with the index fi labeling the diagram. 
Similarly, we obtain for the diagram h 



(2.18) 



T Ilk 



Z GJi$%p(iu] h iCl k ) 



(2.19) 



with 

B ap(iu>i,in k ,<l) = -2T aj (n J 'D(iu) l +in k ,q)np'D(iu l ,q)). 

(2.20) 

The tensor T ai appearing in l|2.20|l describes the renor- 
malization of a current vertex connecting two electronic 
lines with opposite signs of frequencies, 



T af3 = 2 {n a Vn p )\ = 



9 2 

e VpV 



1 + UJ 2 T 2 \ UcT 1 



1 —UJ n T 



(2.21) 



We turn now to diagrams / and g. The expressions for 
the corresponding contributions read 



Sa f a3 (in k ) 



T Q k 



(2.22) 



{fi k + uji)<$> f al3 (-iuji,ini 

-n fe <o)(<o 

Y (Ofc - ui)<£> f a0 (itoi,itl k ) 

W=0 



Sa 9 a (iQ, k ) = — -— — 



4oqT_ 
r fl k 



-n fc <W(<o 



(2.23) 



with 
if 



B T a0 (w h in k ,q)=T a7 (n-yV(iuji+in k ,q)ns)T S p. (2.24) 

The sum of the contributions / and g is therefore given 
by 



p( lUJ h i^lk) 



(2.25) 



We see that when the diagrams / and g are combined, 
the same Matsubara structure as for other diagrams 
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[Eqs. I|2.17|l . I|2.19[l ] arises. In other words, the role of 
the diagrams g is to cancel the extra contribution of di- 
agrams /, which has a different Matsubara structure. 

A word of caution is in order here. In our calculation 
we have set the value of velocity coming from current 
vertices to be equal vf, thus neglecting a particle- hole 
asymmetry. If one goes beyond this approximation and 
takes into account the momentum-dependence of velocity 
(violating the particle-hole symmetry), the above cancel- 
lation ceases to be exact and an additional term with a 
different Matsubara structure arises in Sa f a ^ 9 . After the 
analytical continuation is performed, the corresponding 
correction to the conductivity has a form 



replacement r — > t s 



B lp{iuJi,ittk,ci) = -—T ai T 10 CD(iui,q)V(q,iuJi)), 

ZT S 



and 



(2.29) 



r^inel 
d(J af3 = 



2ao 

T 

d 2 q 



2vr 2Tsinh 2 (w/2T) 



~\ f~ - 7) T ai T lP {V{iu h ci)V{iu u q)). 

(2.30) 

We see that although each of the expressions l|2.28[) . 
(|2.29|) , and (|2.3UI) depends on r s , the single-particle time 
disappears from the total contribution of the Hikami-box, 

a+b+c 1 



-^5BUu,q)lmU{u,q), (2.26) £«+ + c (wfi fe) q) = ^ <2?(t Wj , q)V{iw u q)) 



characteristic for effects governed by inelastic scattering. 
This contribution is determined by real inelastic scatter- 
ing processes with an energy transfer 10 < T and behaves 
(in zero magnetic field) as e 2 (Tr) 2 . This implies that the 
corresponding resistivity correction, Sp ~ (T/eEp) 2 is in- 
dependent of disorder. However, such a correction should 
not exist because of total momentum conservation. In- 
deed, an explicit calculation (see Appendix^) shows that 
this term is canceled by the Aslamazov-Larkin-type dia- 
grams analogous to those describing the Coulomb drag. 

Finally, we consider the diagrams a,b, and c. Already 
taken separately, each of them has the expected Mat- 
subara structure (contrary to the diagrams d,e and f,g, 
which should be combined to get this structure). How- 
ever, another peculiarity should be taken into account. 
The diagrams a, 6, and c form together the Hikami box, 
so that their sum is smaller by a factor ~ t s /t than sepa- 
rate terms. Therefore, some care is required: subleading 
terms of order r s /r should be retained when contribu- 
tions of individual diagrams are calculated. The result 
reads 



4<tq T 

T Q k 



E 

.wj=0 



(2.27) 



Here the contributions of individual diagrams a, b, and 
c have the form 



1. 



f 



5js + (T 1 ) 1 s 



T 5 p{V{iu u q)V(iw u q)) 
(2.28) 

where the matrix T a p has the same form as T a f3 with a 



+ -T ai T 10 {V(iw u q)) , 



(X>(<wi,q))T 7/9 . (2.31) 



The total correction to the conductivity tensor is 
obtained by collecting the contributions (|2.17|) . H2.19f) . 
(|2~2"5l) . and l(2~27jl . Carrying out the analytical continua- 
tion to real frequencies, we get 



Sa a p(Q) 



(TO 



du> lo coth — 
2T 



x [<f> a [3(u + Cl,n)-<f> af3 (u,n)}, (2.32) 



where 



0> Q/3 ( w ,r!) = $«+ b+c (a; ) 0) + ^(a;,fi) 

+ ^(w.fij+^w.n). (2.33) 

We are interested in the case of zero external frequency, 
ft — > 0, when Eq. 1)2. 32[) can be rewritten as 



8 a, 



a/3 



ITTT 



du> 



x $ q/3 (cj,0) — 



1 w ' 
lo coth — 

2T 



(2.34) 



Recalling the definition (|2.18|l of we finally arrive at 
the following result 



8a 



a0 



-2e 2 v 2 F v 



f°° did 


d 








u> coth — 






2TJ 



d 2 q 



Im[ U(u,q) B af} (uj,q) ] , 



(2.35) 



where the tensor B a /s(uj,q} is given by 



B a p(u>,q) = ^(VV) + T aj (^(V)-(n^Vn s )\Tsp 
- 2T ai {n 7 T>npT>) - {T>n a T>n V). (2.36) 
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The first term in l|2.36[l originates from the diagrams 
a,b,c, the second term from a,f,g, the third term from 
h, and the last one - from d and e. We remind the reader 
that this result has been obtained under the assumption 
r s <C t, T _1 ; generalization to arbitrary t s /t and t s T 
will be considered in Sec. Ill Bl It will be shown there 
that the conductivity correction retains the form (|2.35|) 
in the general case but the expression (|2.36[) for B a p{w, q) 
is slightly modified. 



B. General case 

In the previous subsection we have derived the formula 
for the correction to the conductivity tensor for the case 
of a smooth disorder (with r s <C t) assuming r s <C T _1 . 
Since characteristic momenta q and frequencies lo are set 
by the temperature, this assumption implies ql s <C 1 and 
ljt s <C 1. This allowed us to simplify the calculation by 
neglecting the q and u> dependence of Green's functions 
connecting ballistic diffusons and by considering only the 
ladders with many impurity lines. Furthermore, we have 
used the small-angle nature of scattering when calculat- 
ing the Hikami box contribution (|2.31() . We are now go- 
ing to discuss the general case of arbitrary r s /r and Tt s . 

It turns out that the expressions l|2.17l) . I|2.19|) . and 
(|2.25l) for the contribution of the diagrams d — h derived 
in the case of a smooth disorder remain valid in the gen- 
eral situation. The simplest way to show this is to use 
the following technical trick (cf. Refs. 38 39). One can 
add to the system an auxiliary weak smooth random po- 
tential with a long transport scattering time t»t but 
short single-particle f s <C t s , such that Tf s <C 1. This 
potential will not affect the quasiclassical dynamics and 
thus should not change the result. On the other hand, it 
allows us (in view of the condition Tf s <C 1) to perform 
the gradient and frequency expansion in Green's func- 
tions as was done in Sec. Ill Al Adding such an auxiliary 
disorder amounts to a re-distribution between quantum 
and quasiclassical degrees of freedom: all the information 
about the real disorder is now contained in the ballistic 
propagators. It can be verified by a direct calculation 
(without using the additional disorder) that the above 
procedure yields the correct result. 

It remains to consider the Hikami-box contribution 
(|2.27() . When calculating it in Sec. Ill Al we used the 
small- angle nature of scattering implying that a sin- 
gle scattering line inserted between two ballistic prop- 
agators approximately preserves the direction of veloc- 
ity, (VWV) -> {2ttv)- 1 {VV)/t s and (Vn Q Wn p V) -> 
{2-kv)~ 1 (VT>)(1/t s — 1/t)(5 q/ 3. In the more general situ- 
ation, when the scattering is at least partly of the large- 
angle character, this is no longer valid and Eq. (|2.31|) 
acquires a slightly more complicated form, 

B a + b+c {^m,iQ k ) = 
TTvT ai [{VS 1& V) - 2(Vn 7 WnsV)]T Sf 3 



-T ai (V)T~ (l3 , 



(2.37) 



where S xx = S yy = W(n, n'), S xy = -S yx = uj c /2ttu. 

Summarizing the consideration in this subsection, in 
the general situation the interaction correction retains 
the form i|2.35fl with the tensor B a p(uj,q) given by 

B aP (cj,ci) = T ai TTv[{VS lS V) - 2(Vn 1 Wn s V)}Ts0 
-2T ai ( ni V7ipV) - {Vn a Vn p V). (2.38) 



The correction 6p a p to the resistivity tensor is then im- 
mediately obtained by using Sp — —pSa p. This yields 



5p a p = 



9 2 

e VpV 



r°° dej 


d 


LU ' 






& coth — 




doj 


2T\ 



d 2 q 
(2^F 



Im 



(2.39) 

where the tensor B^l is related to B a [3, Eq. Ij2.38|l . via 

#S = (T~ 1 ) ai B 1 s(T~ 1 )8f3. (2.40) 

Explicitly, corrections to the components of the resis- 
tivity tensor are expressed through 5a xx — Sa yy and 
Scr xy = —Sa yx as follows 

Sp xx = pI[{u 2 c t 2 - l)5cr xx + 2w c TS(T X y], (2.41) 
Sp xy = pI[(uj 2 t 2 - l)5a xy - 2(J c t5(t xx \. (2.42) 

Note that the results (|2~^jl . fTSfy for B a p(cj,q) sat- 
isfy the requirement 



B a0 (uj,O) =0, 

as follows from 

(n a T>np)\ Q = a a p{uj)/2e 2 vv 2 F 



(2.43) 



and (2tt)- 1 J dcf> V(<j),(j}') 



9 =o 



i/oj. The condition 



implies that spatially homogeneous fluctuations i n the 
potential do not change the conductivity, see Refs. Il8l42l 
for discussion. 



C. Limiting cases 

Having obtained the general formula, we will now 
demonstrate that it reproduces, in the appropriate lim- 
its, the previously known results for the interaction cor- 
rection. Specifically, in Sec. Ill G"T1 we will consider the 
diff usive li mit Tr < 1 studied in Refs. ITIITolfTTL while 
Sec. lII C 2l is devoted to the B — > case with a white- noise 
disorder addressed in Refs. ll8lT^ . In Sec. Ill C~3l we will 
analyze how the linear-in-T asymptotics of 8<r(B = 0) 
in the ballistic regime obtained in Ref. for a white- 
noise disorder depends on the character of the random 
potential. 
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Diffusive limit 



2. B — * 0, white-noise disorder 



We begin by considering the diffusive limit Tr <C 1 in 
which we reproduce (for arbitrary B and disorder range) 
the logarithmic correction 1)1.20 . (|1.3[) determined by the 
diagrams a-e. Let us briefly outline the corresponding 
calculation. The propagator for ql, lot <C 1 can be de- 
composed as T> = V s + T> r ° s , where V s is singular, while 
2? rcg is finite (regular) at q, to — > 0, see e.g. Refs. liollill 
The singular contribution is governed by the diffusion 
mode and has the form [see Eq. (|U> . 5|> ] 



V s (w,q; <j>,4f) 



V R (<t>, q)^ L W>', q) 



Dq 2 



ILO 



(2.44) 



*„(</>, q) = 1 - ic [ J-> cos(4> - <j) q ) - icl z> sin(</> - <j> q ), 

where D = Vpr/2(l + lo 2 t 2 ) is the diffusion constant in 
the presence of a magnetic field and 



We allow now for arbitrary Tr but consider the limit 
of zero magnetic field assuming a white-noise disorder 
(t = t s and W(n.n') = l/27wr), which is the limit 
studied in Refs. lialTa The contribution l(2~37|) of the 
diagrams a, b, c takes for the white- noise disorder the form 



r>a+b+c _ 1 rp 



(C)(C)^+ Wc <OT) &r4 

T 



T 



6/3- 



(2.49) 



Using now the explicit form of the ballistic propagator for 
the case of white-noise disorder and B — > [Eqs. I|B.4J) . 
(ESJ), i|B~TT)|) . |{B~TT|) . and lfTOB|) ] we recover the results 
for Sa^x and 8p xy obtained in a different way in Refs.fl8l 
and [T^, see Appendix FBI 



'(?) 



4%) 



(9) = 



(?) = 



qvpT 

1 + C^2 T 2 ' 

qv F uj c T 2 

1 + CJ 2 T 2 ' 



(2.45) 
(2.46) 



The leading-order contribution of the diagrams a, 6 
and c (that containing two singular diffusons V s ) is ex- 
actly canceled by the part of the diagrams d and e 
with the structure (2? s n Q 2? reg na2? s ), i.e. with one reg- 
ular part of the propagator inserted between two singu- 
lar diffusons, (V s ) = (Dq 2 — ilo)^ 1 . Indeed, in view of 
(n a T> ies np) — \T a pi the latter contribution reduces to 
— ^(D s } 2 T a f3, while the diagrams a, b and c yield 



\(V 8 ?T ai 



h w c e 7 5 

T 



f sp 



3. B — 0, ballistic limit Tr > 1 

In the ballistic limit Tr 3> 1 and for white-noise dis- 
order the result of Ref. (recovered in Sec. Ill C 21 and 
Appendix[5J) yields a linear-in-T conductivity correction, 
8a = (2vV$e 2 1 'it)Tt for the point-like interaction and 
8a = (e 2 /tt)Tt for the Coulomb interaction. The ques- 
tion we address in this subsection is how this behavior 
depends on the nature of disorder [i.e. on the scattering 
cross-section W(0)]. 

In order to get the Tr ballistic asymptotics, it is suf- 
ficient to keep contributions to (|2.38|) with a minimal 
number of scattering processes. Specifically, the propa- 
gator T> in the first and the third terms of (|2.38(l can be 
replaced by the free propagator, 



(V s YT a0 , 



(2.47) 



where e a( 3 is the antisymmetric tensor, e xx — e yy — 

0) ^xy — ^yx 1- 

It remains thus to calculate only the contribution of 
the diagrams d + e with three singular diffusons, 



2, ,2 



e v 



2- 



oo g 

dui 

— oo 



du> 



lo coth- 



2T 



d 2 q (D s n a V s n V s 



Tm 



(2tt): 
2e 2 v 2 F 



7T(1 + L0 2 c T 2 )- J T 



l + iuo(D s 

l/r 

du> 



d 2 q (-iq a l)(-iQpl) 



(2tt) 2 Dq 2 (Dq 2 - iLof 



2ir 2 



ln(TT)8 a/ 3, 



(2.48) 



V f (iu,q; <j>,<j>') = 



2irS((f> ~ cj)') 



-i(ui + iO) + iqvF cos(<p — <p q ) 



, (2.50) 



while in the second term it should be expanded up to the 
linear term in the scattering cross-section W [the second 
term produces then the same contribution as the first 
term in (|2.38l) ]. The last (fourth) term in 12.38|) does not 
contribute to the Tr asymptotics. We get therefore 

B xx ~ 2nvT 2 ((V i WVi) ~2(V f n x Wn x Vi)) 

- 2r(n x V l n x V i ). (2.51) 

Let us consider first the case of a short-range interaction, 
U (r) = V Q . The structure of Eqs. l(2~351> . (|2~5T|) implies 
that the interaction correction is governed by returns of a 
particle to the original point in a time t < T _1 -C t after 
a single scattering event. It follows that the coefficient 
in front of the linear-in-T term is proportional to the 
backscattering probability W(tt) = w(2kp ), 



in agreement witl 



The result for a point-like in- 



teraction differs only by a factor vVq. 



8 a w 



2vV e 2 



2ttuW(tt)Tt' 2 



(2.52) 
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As shown in Appendix El this result remains valid in 
the case of Coulomb interaction, with the factor 2vVq re- 
placed by unity. This shows that in the ballistic limit 
the Coulomb interaction is effectively reduced to the 
statically screened form, U(r) = l/2v when the lead- 
ing contribution to 5a xx is calculated. According to 
l|2.52|) . in a smooth disorder with a correlation length 
d > k~p} the Tt contribution is suppressed by an expo- 
nentially small factor 2-KVTw(2kp) ~ e~ kpd . In fact, for 
a smooth disorder the linear term represents the lead- 
ing contribution for Tt s 3> 1 only. In the intermedi- 
ate range t _1 <C T <C t^ 1 the dominant return pro- 
cesses are due to many small-angle scattering events. 
However, the corresponding return probability is also 
exponentially suppressed ~ exp(— const r/t) for rele- 
vant (ballistic) times t <C r, yielding a contribution 
5a xx ~ exp[— const (Tt) 1 ' 2 ]. Thus, the interaction cor- 
rection in the ballistic regime is exponentially small at 
B = for the case of smooth disorder. Moreover, the 
same argument applies to the case of a non-zero B, as 
long as 43 uj c < T. 

In any realistic system there will be a finite concen- 
tration of residual impurities located close to the elec- 
tron gas plane and inducing large-angle scattering pro- 
cesses. In other words, a realistic random potential can 
be thought as a superposition of a smooth disorder with 
a transport time r sm and a white-noise disorder charac- 
terized by a time T wn . Neglecting the exponentially small 
contribution of the smooth disorder to the linear term, 
we then find that the ballistic asymptotics (|2.52[1 of the 
interaction correction takes the form 



e t 
5a = Tt x 

7T T wr , 



2i/ Vb, 
1, 



point — like, 
Coulomb, 



(2.53) 



where r _1 — t~^ + t~^ is the total transport scattering 
rate. If the transport is dominated by the smooth dis- 
order, r wn 3> T sm , the coefficient of the Tt term is thus 
strongly reduced as compared to the white-noise result 
ofRef.Ql 

Finally, it is worth mentioning that in addition to the 
Tt term corresponding to the lower limit w ~ T of 
the frequency integral in i|2.35[l . there is a much larger 
but T-independent contribution 5a cx EpT governed by 
the upper limit uj ~ Ep. This contribution is just an 
interaction-induced Fermi-liquid-type renormalization of 
the bare (noninteracting) Drude conductivity. 



III. STRONG B, SMOOTH DISORDER 

A. Quasiclassical dynamics 

We have shown in Sec. Ill CI that due to small-angle 
nature of scattering in a smooth disorder the interaction 
correction is suppressed in the ballistic regime Tt 1 
in zero (or weak) magnetic field. The situation changes 
qualitatively in a strong magnetic field, lu c t 3> 1 and 



lo c 3> T . The particle experiences then within the time 
t ~ T _1 multiple cyclotron returns to the region close 
to the starting point. The corresponding ballistic propa- 
gator satisfies the equation l|2.8ll with the collision term 

The solution of this equation in the limit of a strong 
magnetic field, ui c t 3> 1, is presented in Appendix iDl For 
calculation of the leading order contribution to 5a xx and 
5p xx , the following approximate form is sufficient: 



T>(w, q; (j), 4>') — exp[— iqR c (su\<j) — sin </!/)] 
x(<t>W) e *«W-*') 



E 



Dq 2 — iuo Dq 2 — i(u> — nu> c ) + n 2 /r 

= P s (o;,q;0>')+I? ros (a;,q;0>'), 



(3.1) 



where x(4>) — 1~ iqRc cos 4>/ui c t andD ~ R 2 /2t, and the 
polar angles of velocities are counted from the angle of q. 
Equation 1)3. lfl is valid under the assumption (qR c ) 2 <C 
lu c t. We will see below that the characteristic momenta 
q are determined by the condition Dq 2 ~ lo ~ T, so 
that the above assumption is justified in view of u c T. 
Furthermore, this condition allows us to keep only the 
first (singular) term D s in square brackets in l|3.1|l when 
calculating (T>), 



(V) 



Dq 2 - i 



(3.2) 



where Jq(z) is the Bessel function. Moreover, the formula 
(|2.36(l for B xx can be cast in a form linear in T> by using 



(£>£>) 
(n a VnpV) 

(T>n x 'Dn x 'D) 



i d 
vf dqp 
1 d 2 



2v% 



n a V), 



(3.3) 
(3.4) 

(3.5) 



Therefore, it is again sufficient to take into account only 
the first term in 1)3. ljl for calculation of B xx if the identi- 
ties (|3.3ll . H3.4fl . and (|3.5|l are used. (Of course, B xx can 
also be evaluated directly from Eq. I|2.36|l . but then the 
second (regular) term T) ICS in l|3.1|l has to be included.) 
Combining all four terms in (|2.36(l . we get 



J 2 (gR c ) D T q 2 
(co c t) 2 (Dq 2 -iujf 
4r 3 J 2 (Q)Q 2 
(3 2 (Q 2 -iQ) 3 ' 



(3.6) 



In the second line we introduced dimensionless variables 
Q = qR cl f2 = 2lut, (3 = lu c t. 

Note that Eqs. I|3.2() . (|3.6I) differ from those obtained 
in the diffusive regime by the factor J^(qR c ) only. This is 
related to the fact that the motion of the guiding center is 
diffusive even on the ballistic time scale t <C r (provided 
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t 3> w^ 1 ), while the additional factor corresponds to the 
averaging over the cyclotron orbit (see Sec. UVI below). 

We turn now to the calculation of B xy . Substituting 
(|3.1|) in (|2.36|) , we classify the obtained contributions ac- 
cording to powers of the small parameter 1//3. The lead- 
ing contributions are generated by the first and the last 
terms in H2.36|l and are of order 1/(3, i.e. larger by factor 
(3 as compared to B xx , Eq. (|3.6p . (This extra factor of 
(3 is simply related to \o xy \fa xx — (3.) However, these 
leading contributions cancel, 



^(VV) - (Vn x Vn y V) 



order 1//3 



— (V S V S ) - {V s n x V^n y V s ) 
2r 3 J 2 (Q) , 2r 3 J 2 (Q) 



(3.7) 



p {Q 2 -iny (3 {Q 2 -iny 



o, 



as in the diffusive limit, see the text above Eq. I|2.47|l . 

To evaluate terms of higher order in 1//3, we need a 
more accurate form of the propagator l|3.1|) . Since the 
contributions of order 1//3 2 to B xy turn out to cancel 
as well, we have to know the propagator with the accu- 
racy allowing to evaluate the terms of order 1//3 3 . To 
simplify the calculation, we use again the identities 1)3. 3|1 
and 1)3. 4JI . As to Eq. (|3.5(l , it cannot be generalized onto 
xy-componcnt of the tensor, and we use instead 



% I dV 
{Vn x Vn y V) — — ( -^—n y V 



(3.8) 



It is then sufficient to calculate the propagator T> up 
to the 1//3 2 order. This is done in Appendix [D] see 
Eqs. (|D.14|) - (|D.17|I . Substituting this result for V in 
Eq. H3.6fl an d combining all terms, we get after some 
algebra 



B xy (uj,q) = --^ 



7Q 2 J 2 (Q)+4QJ (Q)Ji(Q) 



(O 2 - iny 



4QMQ)MQ)) 

Q 2 -ifl 



(3.9) 



We see that similarly to l|3.6[) the kernel B xy (uj 7 q) has 
a diffusive-type structure with Q 2 — ifl in denominator 
reflecting the diffusion of the guiding center, while the 
Bessel functions describe the averaging over the cyclotron 
orbit. Clearly, both kernels l|3.6|l and (|3.9|l vanish at 
q = 0, as required by 



B. Point-like interaction 

To find the interaction correction to the conductivity, 
we have to substitute Eqs. I|3.6I) and l|3.9|) in the formula 
(|2.35|l . We consider first the simplest situation, when 
the interaction U(oj, q) in (|2.35|) is of point-like character, 
U(uj, q) — Vq. Using Vpqdq = u} 2 QdQ, we see that all the 



B-dependence drops out from 5a xx , and the exchange 
contribution reads 



Sa x 



-8e 2 vV Q 



dtt d 



fl coth- 



"47V 



QdQ J 2 (Q)Q 2 

1m- 

2tt 



(Q 2 -in)f 



(3.10) 



To simplify the result (|3.10|) . it is convenient to perform 
a Fourier transformation with respect to £1 (which corre- 
sponds to switching to the time representation) 



Im 



— F(w) — 
nT 2 t 



uj coth — 
2T 



dt 



F(t). 



lo sinh 2 (7rTt) 
The integral over Q is then easily evaluated, yielding 



(3.11) 



S<T XX 

G (x) 



■—uV G (Tt), 



2 2 
7T X 



' dueicp(— 1/u) 



'o sinh 2 (ttxu) 
x [(«-l)7 (l/u) + Ji(l/«)], 



(3.12) 
(3.13) 



where Iq{z) and h{z) are modified Bessel functions. The 
Hartree term in this case is of the opposite sign and twice 
larger due to the spin summation (we neglect here the 
Zeeman splitting and will return to it later). 

It follows from Eqs. I|3.6|l and l|3.9|l that the correction 
to the Hall conductivity is smaller by the factor (cj c t) _1 
as compared to (|3.12H . This implies, according to (|2.41|) 
that in a strong magnetic field the correction to the lon- 
gitudinal resistivity is governed by Sa xx , 



Spx 



pa 



(TO 



(3.14) 



similarly to the diffusive limit (|1.3|) . In fact, it turns out 
that the relation (|3.14|l holds in a strong magnetic field, 
lj c 3> T, for arbitrary disorder and interaction, see below. 
On the other hand, as is seen from (|2.42|) . contributions 
of both Sa xx and 5a xy to Sp xy are of the same order in 
(lUct)^ 1 . We will return to the calculation of Sp xy in 

secnnm 

The MR p xx (B) is thus quadratic in B, with the tem- 
perature dependence determined by the function Go (TV) , 
which is shown in Fig.0^,. In the diffusive (x <C 1) and 
ballistic (x ^> 1) limits the function Go(x) has the fol- 
lowing asymptotics 



Go (a;) 



with 



c x 



CO 



lax + const, 
- 1 ' 2 , z»l, 



3C(3/2) 
160F 



0.276 



(3.15) 



(3.16) 



11 



o 
0.01 



V 

\ 
\ 

V\ 

\\ 

\\ 

\\ 


— G (x) 


— -1.7-lnx 


\\ 
\\ 
\\ 
» \\ 
\\ \ 
^ \ 
v 4 \ 

V \ 
\\ \ 


-1/? 

— 0.276 x 


\ ~\ 

a \ x 

\ 

i \ i 





teraction 



0.1 1 

Tt 



10 




FIG. 4: Functions G (Tt) (a) and G F (Tr) (b) determin- 
ing the T-dependence of the exchange term for point-like, 
Eq. I|M.12|I . and Coulomb, Eq. Ij3.19ft . interaction, respectively. 
Diffusive and ballistic asymptotics, Eq. (13. 151 and Eq. 13. 211 . 
are also shown. 



(here ((z) is the Riemann zeta- function). Let us note 
that the crossover between the two limits takes place at 
numerically small values Tt ~ 0.1 (a similar observation 
was made in Refs. ll8lll9|h This can be traced back to the 
fact that the natural dimensionless variable in l|3.12[l is 
2-kTt. 



C. Coulomb interaction, exchange 

For the case of the Coulomb interaction the result 
turns out to be qualitatively similar. Substituting l|3.2l) 
in l|2.14[) and neglecting the first term q ~ (T/D) 1 /' 2 <C n 
in the denominator of i|2.14[l . we obtain the effective in- 



o 2 - in 



(3.17) 



2v Q 2 -ifl[l-J 2 {Q)Y 

Inserting l|3.17[) and 1)3. 6fl into (|2.35|) , we get the exchange 
(Fock) contribution 



Sat 



Im 



QdQ 



d_ 

on 



D, coth- 



4Tr_ 
Q 2 J 2 (Q) 



(Q 2 - 



(3.18) 



Using lj3~T4"jl wc find the MR 

Sp F xx (B) (cu c t) 2 



Pa 



■nkpl 



G f {Tt) 



poo 

G F {x) = 32n 2 x 2 / dQQ 3 J%(Q) 
Jo 

^ n (12TTxn[l - J 2 (Q)} + [3 - J 2 (Q)]Q ; 



(3.19) 
(3.20) 



^ {Airxn + Q 2 )3(Airxn[l - J 2 {Q)] + Q 2 ) 2 " 

Note that in contrast to the case of a point-like inter- 
action, a transformation to the time representation does 
not allow us to simplify l|3.18|l . since the resulting Q- 
integral can not be evaluated analytically. We have cho- 
sen therefore to perform the O-integration, which results 
in an infinite sum (|3.20|) . This amounts to returning to 
the Matsubara (imaginary frequency) representation and 
is convenient for the purpose of numerical evaluation of 
Gp(x). In the diffusive (x <C 1) and ballistic (x ^> 1) 
limits this function has the asymptotics 



G F {x) ~ 
and is shown in Fig. 0Jd. 



— In a; + const, 

£° -1/2 

2 X ' 



X < 1, 
X > 1, 



(3.21) 



D. Coulomb interaction, Hartree contribution 

We turn now to the Hartree contribution. The corre- 
sponding diagrams can be generated in a way similar to 
exchange diagrams (Sec. Ill All but in this case one should 
start from two electron bubbles connected by an interac- 
tion line. There are again two distinct ways to generate a 
skeleton diagram: two current vertices can be inserted ei- 
ther both in the same bubble or in two different bubbles. 
Then one puts signs of Matsubara frequencies in all pos- 
sible ways and insert ballistic diffusons correspondingly. 
The obtained set of diagrams is shown in Fig. [3] There 
is one-to-one correspondence between these Hartree di- 
agrams and the exchange diagrams of Fig. ^ which is 
reflected in the labeling of diagrams. 

As seen from comparison of Figs. and [5] the elec- 
tronic part S„«(<^, 4>') of each Hartree diagram is identi- 
cal to that of its exchange counterpart. The only differ- 
ence is in the arguments of the interaction propagator, 
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FIG. 5: Hartree diagrams for the interaction correction to 
a a 0. The diagrams are labeled in the way as their exchange 
counterparts in Fig. The diagrams obtained by a flip 
and/or by an exchange + <-> — should also be included. 



U(ui,q) — * U(0, 2fcjrsin ), where and <j>' are polar 
angles of the initial and final velocities [cf. Eqs. I|2.11[l . 
(|2.12|) ]. Therefore, in the first order in the interaction, 
the Hartree correction to conductivity has a form very 
similar to the exchange correction (|2.35() . 



So^p = 4e 2 v%v 





d 








ui coth — 


Loo 2^ 


duj 


2TJ 



d 2 q [ d(f> #' 
(2tt) 2 J 2^2^ 
Im [Un^^Ba^^^A') 



(3.22) 



where 



U 0,2fc F sin- 



(3.23) 



is the Hartree interaction and B a p(v, q; 0, <fi') is given by 
Eqs. (|2.3(il) . (|2.38l) without angular brackets (denoting 
integration over and </>'), see Eq. (I2.12|) . Clearly, for a 
point-like interaction U(ui, q) = Vq this yields 



(3.24) 



as has already been mentioned in Sec. IIII Bl 

In the case of the Coulomb interaction the situation 
is, however, more delicate 4 ^. To analyze this case, it is 
convenient to split the interaction into the singlet and 
triplet part o 1 ! 18 ! 44 . For the weak interaction, k -C hp, 
the conductivity correction in the triplet channel is then 
given by Eq. I|3.22|l with an extra factor |. 

As to the singlet part, it is renormalized by mixing with 
the exchange term. The effective interaction U s in the 
singlet channel is therefore determined by the equation 

Z J Zn Ztt 



x [E/o - ^(01,02)C/ S (02,0'), 

(3.25) 

where Uq = 2ire 2 /q is the bare Coulomb interaction, and 

V(uj, q; 4>i, 2 ) = 2v[2t:5{4>i - 2 ) + iujT>(uj, q; 01,02)] 

(3.26) 

describes the electronic bubble. Solving (|3.25(1 to the first 
order in £/h, we get 

U\w, q; 0, 0') = U{u>, q) - U%(w, q; 0, 0'), (3.27) 

where t/(u;,q) is the RPA-screened Coulomb interaction 
(|2.14ll which has already been considered in Sec. IIII CI 
while the second term describes the renormalized Hartree 
interaction in the singlet channel, 



,/s 1_ f d<pi d(p 2 

' 2n / 2tt 2tt 



X [/J H (0, 0l)P(01, 02) + P(01, 2 )t/ H (02, 0')] 

1 f d<fii d4>2 d<p3 dcj>4 

+ 2ip J ^T^T^r^r 

X P(0i,0 2 )f/ H (02, 03)^(03, 04). 



(3.28) 



Here n = (P) is the polarization operator (|2.15|) . and 
we have used the singular nature of the bare Coulomb 
interaction implying |n|C/o S> 1 for all relevant momenta. 

Taking into account that the angular dependence of 
leading contributions to B xx (u>, q; 0, 0') and T>(u>, q; 0, 0') 
is of the form exp[— iQ(sin0 — sin0')], we find that the 
singlet Hartree correction to a xx is given by Eq. (|3.22|l 
with a replacement 



(^(0,00) I w 

4[l + iw(P(w,q))] 2 



(3.29) 



Note that in the diffusive limit B a p is independent of 
0, 0', so that only the zero angular harmonic of the inter- 
action contributes. On the other hand, the zero angular 
harmonic is suppressed in the effective singlet- channel 
interaction (|3.29|l . Therefore, the singlet channel does 
not contribute to the Hartree correction in the diffu- 
sive limit, in agreement with Refs. IlKil The situation 
changes, however, in the ballistic regime, when B a p be- 
comes angle- dependent. 

After the angular integration, the triplet Hartree con- 
ductivity correction takes the form 13.10fl with the re- 
placement Vq — > l/2u, and 



Jo(Q) 



-3y 



dtp Jo(2Q sin< 
27r y + 2 sine 



(3.30) 



where y = n/kp. For the singlet part we have a result 
similar to (|3.18() with a slightly different Q-integral, 



QdQ 



J(y,Q)Q 2 



(Q 2 -in[i~ j 2 (Q)]) 2 (Q 2 -ir!)' 
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where 



J(y,Q) = -,r? M2Qsia !\- JS(Q) . (3.31) 



/q 2tt y + 2 sin 
This yields for the total Hartree contribution 



■[G s li (Tr,y) + 3G\ i (Tr,y% (3.32) 



where Gfj and G^ governing the temperature dependence 
of the singlet and triplet contributions have the form 

/>oo 

G s n (x,y) = 32ir 2 x 2 / dQQ 3 J(y,Q) 
Jo 

^ n(12irxn[l - J 2 (Q)] + [3 - 2J 2 (Q)]Q 2 ) 
X hi ( i7:xn + Q 2 ) 2 (^xn[l - J 2 (Q)] + Q 2 )3 ' 

(3.33) 

7TX 7/ 



4 J sinh 2 (-row) 

* n exp[-2sin 2 0/u] 

dtp ■ 

y + 2 sin 



(u — 2 sin 2 



(3.34) 



Figure EJi shows Gh(x,J/) = G|j(^,y) + 3G^(a;,j/) as a 
function of x = Tt for several values of y = n/kF- The 
asymptotic behavior of Sp^ x is as follows: 



5p« (B) (co c r) 2 



nkpl 

' y^-y 



ln(Tr) - 
tflr^yCZY) 1 /*], 
[ ircoiTr)-^, 



lny], Tr«l, 

1 < Tt < 1/y 2 , 

Tr » 1/y 2 , 
(3.35) 



We see that at ft/fc^ 1 a new energy scale Th ~ 
T~ 1 (kp / k) 2 arises where the MR changes sign. Specif- 
ically, at T « T H the MR, Sp xx = Spl x + Sp^ x , is 
dominated by the exchange term and is therefore neg- 
ative, while at T » Tjj the interaction becomes ef- 
fectively point-like and the Hartree term wins, 5p^ x = 
—2Sp xx , leading to a positive MR with the same (Tr) -1 / 2 
temperature-dependence, see Fig.EK- 



E. Hartree contribution for a strong interaction 



in sec. rnrni we have assumed that n/hi? <C 1, or, in 
other words, the interaction parameter r s = \J~2e 2 /evf 
(where e is the static dielectric constant of the mate- 
rial) is small. This condition is, however, typically not 
met in experiments on semiconductor structures. If k/Izf 
is not small, the exchange contribution (|3.19|) remains 
unchanged, while the Hartree term is subject to strong 
Fermi-liquid renormalizationii^ and is determined by 
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FIG. 6: Hartree contribution, Gh{Tt), for (a) weak interac- 
tion, n/kF = 0.1, 0.2, 0.3, 0.5, and (b) strong interaction, 
Fo = —0.3, —0.4, —0.5 (from bottom to top). Dashed curves 
represent the exchange contribution. 



angular harmonics F^ p of the Fermi-liquid interaction 
F a ' p {0) in the triplet (a) and singlet (p) channels. 

The effective interaction U"^ replacing [/h(</>, ^') in 
(|3.22() is then given by an equation of the type (|3.25() 
but with F a ' p ((t) - (f>')/p substituted for U H (^,0') (and 
without Uq in the triplet channel), 



E>2f(0.0') 



FP{ 



2v 



u? s (<t>,<f>') = 



C/ + 
Uq 

F a {4> — 4>') 



~2ttYk 



2v 



(3.36) 



2v 



(tyi d<p 2 
~2n~2n 



F a ( 



2v 



T(<p l7 cb2)U° s (fo,<b'), (3.37) 



A general solution of these equations requires inversion 
of integral operators with the kernels / — F a V and I — 



14 



(Uq + F p yP and is of little use for practical purposes. 
The situation simplifies, however, in both diffusive and 
ballistic limits. 

In the diffusive regime, T <C 1/t, the second term in 
the polarization bubble l|3.26l) and B a p are independent 
of angles 4», <p' . As discussed in Sec. 1111 Dl this leads to 
the suppression of the Hartree contribution in the singlet 
channel, while in the triplet channel only the zero angular 
harmonic survives, 



1 F£(Dq 2 -iuj) 
2^ (1 + F°)Dq 2 - iuj 



(3.38) 



We then reproduce the known resurfci*— Gh(Tt) 
3G^(Tr) with 



£=0 
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GUTt) 



In Tt. 



(3.39) 



In the ballistic limit, T 1/t, the first term is 
dominant in (|3.26ll . since (T>) is suppressed by a factor 
Jq(Q) -C 1, according to l|3.2l) . The angular harmon- 
ics then simply decouple in Eqs. (|3.3tj[l . I|3.37|l . yield- 
ing effective Hartree interaction constants E/g cff = 0, 
K, e s = (2^)- 1 ^/(l + m ± 0, and U^ eS = 

(2vy 1 F^ n /{l + F^ n ). Therefore, the Hartree contribution 
reads 



GhP>) = -| 



m/0 



FP 



3 St 



per 



po 



1 



(3.40) 

From a practical point of view, it is rather inconvenient 
to describe the interaction by an infinite set of unknown 
parameters F% p . For this reason, one often assumes that 
the interaction is isotropic and thus characterized by two 
coupling constants Ffi and Ffi only. Within this fre- 
quently used (though parametrically uncontrolled) ap- 
proximation, the singlet part of the Hartree term is com- 
pletely suppressed. The Hartree contribution is then de- 
termined solely by the triplet channel with the effective 
interaction 



1 F n ° 



Q 2 - iil 



2v 1 



per 

Q 2 



i 



1 + Ft 



:Jo(Q) 



The Hartree correction to the resistivity takes the form 
of Eq. (|3.19|l with an additional overall factor of 3 and 
with J 2 (Q) multiplied by F CT /(1 + F CT ), 



Jo(Q) -> Jo(Q) 



Fn 



(3.42) 



everywhere in (|3.20|) ; the result is shown in Fig. HjJ} for 



several values of FZ . 



FIG. 7: The function G^{TT,e z ), Eq. iTHoTl . describing the 
temperature dependence of the triplet contribution is shown 
for Fq — —0.3 and different values of Zeeman splitting, 



e z = 0.1, 0.3, 0.5, 1.0 (from top to bottom), 
represents the case e z = 0. 



Dashed curve 



F. Effect of Zeeman splitting 

Until now we assumed that the temperature is much 
larger than the Zeeman splitting Ez , T ^S> Ez ■ In typical 
semiconductor structures this condition is usually met in 
non-quantizing magnetic fields in the ballistic range of 
temperatures, allowing one to neglect the Zeeman term. 
If, however, this condition is violated, T < Ez, the Zee- 
man splitting suppresses the triplet contributions with 
the z-projection of the total spin S z = ±1, while the 
triplet with S z = and singlet parts remain unchanged. 

In the case of a weak interaction, n/kp <C 1, the triplet 
contribution 3G#(Tt, n/kp) in Ea. 13.32fl is modified in 
the following manner, 

3G^(x, y) -> G^x, y) + 2Re & K (x, y; e z ), (3.43) 

where e z = 2tEz, and the function G^(a;,?/;e z ) describ- 
ing the temperature dependence of the contribution with 
±1 projection of the total spin is given by 



G^(a;,y;e z ) 



irx 2 y 



' du exp[ie z u] 
o sinh 2 (7ra;it) 
exp[— 2 sin 2 tf>/u] 
y + 2 sin <p 



2 sin 2 



(3.44) 



We see that at Tt <C e z , the contributions of ±1- 
components of the triplet saturate at the value given by 
(|3.34|) with a replacement Tt — > e z , i.e. at ~ G^(e z ,y). 
In the opposite limit, Tt 3> e z , we have G|j(x, y;e z ) ~ 
Gjj(x, 2/), and the result l|3.32|l is restored. 

The triplet contribution for strong isotropic interaction 
(i.e. determined by Fq only) in the presence of Zeeman 
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splitting reads 
5pl x {B) _ Kt) 2 



Po 



■nkpl 



[GUTr,0) + 2Re GUTr,e z )}. 



(3.45) 



The function Gy(Tr, e z ) is given by a formula similar to 

6129, 

roc 

G£(x, e z ) = 32ir 2 x 2 / dQQ 3 [l - J a {Q)) 
Jo 

n(12TrxnJ°(Q) + [2 + J'(Q)][Q 2 + ie z \) 



J (4^xn + [Q 2 + ie z ]) 3 (47ra;nJ CT (Q) + [Q 2 + ie,]) 2 ' 

(3.46) 

with J°(Q) as defined in ((3.421) . Again, for high temper- 
atures Tt 3> e z , all the triplet components contribute, so 
that the overall factor of 3 (as in the absence of the Zee- 
man splitting) restores. On the other hand, for Tt -C e z , 
the contributions with ±1 projection of the spin saturate 
at low temperatures, and therefore the triplet contribu- 
tion is partly suppressed, see Fig. [7| 



G. Hall resistivity 



As discussed in Sec. lIIIBl calculation of the correction 
8p xy to the Hall resistivity requires evaluation of both 
So xx and Sa xy . In fact, as we show below, the temper- 
ature dependence of Sp xy in a strong magnetic field is 
governed by Sa xx in the diffusive limit and by Sa xy in 
the ballistic limit. 

Since 5a xx has been studied above, it remains to cal- 
culate 5a xy . Using the result ((3.9(1 for the corresponding- 
kernel B xy , we get the exchange contribution for the case 
of a point-like interaction 



e 2 V V Q 

OOxy = -7T^ 

2lT Z U} r T 



G (xv) {Tt)+G 



UV 



(3.47) 



where the temperature dependence of the correction is 
governed by the function 



G xy) (x) = -tt 2 f° — exp(-l/u) 
Jo u 



(3.48) 



_ sinh 2 (7rxu) (ttu) 2 
x [ (9u - 3)I (l/iO + (3 - 2«)7i(l/u) ] . 

When writing H3.48fl . we subtracted a temperature inde- 
pendent but ultraviolet-divergent (i.e. determined by the 

upper limit in frequency integral) contribution G^y ; we 
will return to it in the end of this subsection. 

The function G Q xv \x) has the following asymptotics: 



G^\x) 



9nx, 

llciz 1 / 2 , 



x -C 1, 
x > 1, 



(3.49) 



with 



ci = -\C{l/2) ~ 0.647. 



(3.50) 



Combining 1(3.12(1 and 1(3.47(1 and using (|2.42|l . we find 
the correction to the Hall resistivity, 



Pxy nkpl 



(3.51) 



where 



GS"(x) = 2G (x) - G {xy \x) 

-2 Ins + const, i<1, 



Tlcix 



1/2 



X > 1. 



(3.52) 



The function Gq xb (x) is shown in Fig. |H1 As usual, the 
Hartree term in the case of point-like interaction has an 
opposite sign and is twice larger in magnitude, if the 
Zeeman splitting can be neglected. 

An analogous consideration for the Coulomb interac- 
tion yields a similar result for the exchange correction 



S Ply 
Pxy 

G p ^{x) 



GPx 
F 



(TV) 



2G F (x)-G {xy) (x) 
( —2 In x + const, x <C 1, 



11 1/2 

-ycix , 



x > 1. 



(3.53) 



(3.54) 



The function Gp v (x) is obtained by substituting (|3.17() 
and 1)3.9(1 in 1(2.35(1 [cf. similar calculation for 8cr xx leading 
to Eqs. ipn%j> and (|3~2U|) .] The function G^ v (x) describ- 
ing the temperature dependence of the exchange correc- 
tion to the Hall resistivity is shown in Fig. [H] In the bal- 
listic regime, where G F xy \x) dominates, the interaction 



becomes effectively point-like with i/Vq — i 
can simplify the calculation using G F ' (x) 



so that one 
■\G[ xy \x). 



To analyze the Hartree contribution, we restrict our- 
selves to the isotropic-interaction approximation. Then, 
similarly to the consideration in the end of Sec. IIII El 
only the triplet part contributes, and, in order to cal- 
culate G {xv \x), one should use Eqs. ((3.9(1 and 1(3.41(1 . 
In the diffusive limit the Hartree correction to the Hall 
resistivity is determined by ((3.39|) . while in the ballis- 
tic limit we have again effectively point-like interaction 
with uV = §F °7(1 + F CT ), implying that G^ v) (x) ~ 
-3G t xv \x)F£/2(l + F$). This yields 



S Ply 



GS»(Tt) 



nkfl 



G^ y {x) ~ 3x 



1 - 



ln(l + F%) 



FR 



In. 
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F; 



— Ci 

2 1 + F CT 



v—x 1 ' 2 , 



(3.55) 
x < 1, 

x > 1. 
(3.56) 
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FIG. 8: Functions G P xy (Tr) (lower curve)and Gp y (Tr) 
(upper curve) describing the temperature dependence of the 
Hall resistivity for point-like and Coulomb interaction, respec- 
tively. Diffusive (x <C 1) and ballistic (a; ^> 1) asymptotics, 
Eqs. 113.5211 and (13.5411 . are also shown. 



We return now to the T-independent contribution 
cjjy^ that was subtracted in Eq. H3.48[l . In view of the 



divergency of this term at u 



short-time cut-off u n 



-> 0, it is determined by the 
,/2r 



G 



uv 



du 



-1/2 



(3.57) 



Since the correction we are discussing is governed by cy- 
clotron returns, the cut-off t m i n corresponds to a single 
cyclotron revolution, u m i n ~ tt/ui c t. [On a more formal 
level, this is related to the assumption oj <C oj c used for 
derivation of (|3.48l) : see the text below Eq. (|3.1|) .] We 



have, therefore, G 



_ Jxy) 



uv 

AS 



(wcr) 1 / 2 , with a constant 
of order unity^. For the point-like interaction, 
the considered term produces a temperature-independent 
correction to the Hall resistivity of the form 



5 Pxy = vVo&ri 

Pxy nkpl 



(UJ C T 



,1/2 



(3.58) 



In the case of Coulomb interaction, this correction (with 
both, exchange and Hartree, terms included) has the 
same form with uV \ [1 - 3F CT /(1 + Fg)] . 

Finally, let us discuss the expected experimental man- 
ifestation of the results of this subsection. Equations 
1)3.54(1 . (|3.55() predict that in the presence of interaction 
the temperature-dependent part of the Hall resistivity 
p xy (B) in a strong magnetic field lu c 3> t _1 ,T is linear 
in B at arbitrary T, with the T-dependence crossing over 
from InT in the diffusive regime to T 1 / 2 in the ballistic 
regime. More specifically, if the interaction is not too 
strong, the exchange contribution (|3.54|) wins and the 
slope decreases with increasing temperature, while in the 



limit of strong interaction the slope increases due to the 
Hartree term (|3.55|l . In an intermediate range of Fq the 
slope is a non-monotonous function of temperature. Sur- 
prisingly, this behavior of the slope of the Hall resistivity 
is similar to the behavior of a xx obtained in Ref. [18| for 
B = and white- noise disorder. This is a very non-trivial 
similarity, since the correction to p xy at weak fields^ 
shows a completely different behavior, vanishing as T _1 
in the ballistic regime. In addition to the temperature- 
dependent linear-in-i? contribution, the interaction gives 
rise to a T-independent correction (|3.58l) , which scales as 
5p xy oc B 3 / 2 (assuming again that lu c ^> r _1 , T) 

Let us recall that these results are governed by mul- 
tiple cyclotron returns and thus are valid under the as- 
sumption u c 3> T. In the opposite case, lu c <C T, the 
correction is suppressed in the ballistic regime (similarly 
to 5p xx , see Sees. Ill BT> and llll A)l . and the Hall resistance 
takes its Drude value. 



IV. QUALITATIVE INTERPRETATION: 
RELATION TO RETURN PROBABILITY 

It was argued in Ref. by using the Gutzwiller trace 
formula and Hartree-Fock approximation that the inter- 
action correction to conductivity is related to a classical 
return probability. The aim of this section is to demon- 
strate how this relation follows from the explicit formulas 
for a xx . 

We begin by considering the case of smooth disorder, 
when the kernel B xx (u,q) is given by Eq. I|2.36|l . For 
simplicity, we will further assume a point-like interac- 
tion, when only the first two terms in H'2.36[) give non- 
zero contributions. In fact, we know that the result for 
the Coulomb interaction is qualitatively the same [cf. 
Eqs. (|3~T5fl and EHJ]. 

We will concentrate on the first term in (|2.36|l : the 
second one yields a contribution of the same order in the 
ballistic regime and is negligible in the diffusive limit. 
Therefore, for the purpose of a qualitative discussion it 
is sufficient to consider the first term. Using (|3.3|) . the 
corresponding contribution can be estimated as 



8a ^ 



To 



— oc 
o 

dt 



d 



lu coth- 



du> 
(nT)H 



(dq)Re 



io sinh (nTt) 

V [ dt(D(r = Q,t)) 
Jo 



2TJ 

t(V(r = 0,t)) 



dui 



(4.1) 



where a xx is the Drude conductivity in magnetic field and 
we performed in the second line the Fourier transforma- 
tion of T> to the coordinate-time representation (|3.11l) . 
The return probability in a strong magnetic field, 

LO c T > 1, 



R(t) = (V(r = 0,t)), 



(4.2) 
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Diffusive 




FIG. 9: Schematic plot of the return probability R(t) in 
a strong magnetic field and smooth disorder. In the ballis- 
tic regime, the peaks are separated by the cyclotron period, 
ii = 2n/ui c . Dashed curve represents the smoothened return 
probability R(t) 



is shown schematically in Fig. [5J In the diffusive time 
range, t 3> r, it is given by R(t) — l/4wDt (where D is 
the diffusion constant in the magnetic field, D ~ R^/2t). 
Equation l|4.1|l thus yields in the diffusive regime, Tt <C 
1, 



8&xx Vo 



Or 



D 



ln(Tr), 



(4.3) 



in agreement with l|1.2l) , 12.48(1 . 

At short (ballistic) time, t <C r, the return prob- 
ability is governed by multiple cyclotron returns after 
n = 1, 2, . . . revolutions, 



R W = A K CT > D2 6XP 



[t - 2irn/LU c ] 2 uj3i 
127m 



(4.4) 

Since T <C u> c , the conductivity correction (|4.1|l is in fact 
determined by the smoothened return probability, 



R(t) 



1 1 /r\i/2 



(2tt) 3 / 2 R 2 c 



(4.5) 



Substituting l|4.5|l in i|4.1[l we find that in the ballistic 
limit, Tt ^> 1, the conductivity correction scales as 



D 



(Tt) 



-1/2 



(4.6) 



in agreement with the exact results l|3.12[) . (|3.15() . As to 
the diffusive regime, Tt <C 1, the contribution of short 
times t < t to the integrand in (|4.1|l yields a subleading 
T- independent correction ~ Vo/D to l|4.3|l . 

It is worth emphasizing that the ballistic behavior 
(I4.5|l of the return probability R(t) corresponds to a one- 
dimensional diffusion. Consequently the ballistic result 
(|4.(i|) has the same form as the diffusive Altshuler-Aronov 
correction in the quasi-one-dimensional geometry. To 
clarify the reason for emergence of the one-dimensional 
diffusion, we illustrate the dynamics of a particle subject 
to a strong magnetic field and smooth disorder in Fig. I1UI 




X 



X 



FIG. 10: Schematic illustration of the ballistic dynamics in 
a strong magnetic field. The thick line shows the particle 
trajectory (two cyclotron revolutions disturbed by the smooth 
random potential) . The thin line is the diffusive trajectory of 
the guiding center. 



Let us assume that the velocity is in y direction at 
t = 0. As is clear from Fig. ^3 the return probability 
Ri after the first cyclotron revolution (the integral of the 
first peak in Fig. |5J) is determined by the shift S x of the 
guiding center in the cyclotron period t\ — 2tt/uj c . In 
view of the diffusive dynamics of the guiding center, this 
shift has a Gaussian distribution with 



St 



(81) = 2Dh 



2nR 2 r 



(4.7) 



yielding 



Ri 



(0J c t) 



1/2 



2-kvfRc 



nS 2 



(4.8) 

after n revolutions, 
= Ri/y/n. As to the 



Furthermore, we have (S 2 ) 
yielding the return probability R 
y-component S y of the guiding center shift, it only gov- 
erns the width of the peaks in Eq. (|4.4|l and Fig.0witliout 
affecting R(t). Therefore, the smoothened return proba- 
bility is 



R(t) = -f- , (4.9) 

r l n =t/t! 

which reproduces Eq 
As mentioned in Sec 



the emergence of the one- 
dimensional diffusion in the ballistic regime is reflected 
by the factor Jq(Q) ~ 1/ttQ in the formula l(3.6|l for 
the kernel B xx (oj,q). This factor effectively reduces the 
dimensionality of the g-integral, J d 2 q — ► R~ x J dq. 

In the above we considered a system with smooth dis- 
order, for which 6a xx at B — vanishes exponentially in 
the ballistic limit. Now we turn to the opposite case of 
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a white-noise disorder. We will show that the linear-in- 
T corrections 7 ^ (Sec. Ill CB ) is again related to the re- 
turn probability but the relation is different from 114. 1|) . 
Indeed, according to l|2.49|) . we have now the structure 
{T>)(T>) instead of (TJD) that was relevant for smooth 
disorder. On the other hand, the return probability at 
ballistic times t <C t is clearly dominated by processes 
with a single back-scattering event, implying 



(V(r = 0,f» ~ - / (dq)du(V i (io,c i )) 2 e 



2iut 



(4.10) 



Therefore, the contribution of the first term in 12.49fl can 
be cast in the form 

^ ~ Vor r 1 f(P(r = 0,t)) ) 

- t/orlconst-^^^-T- 1 ))]. (4.11) 

It is easy to see that the probability of a ballistic return 
after a single scattering event is 



(X>(r = 0,i)) 



1 f, 2 8{t-2r 1 /v F ) 1 



(vfti) 



vtrt 



(4.12) 



Substituting (|4.12|) in 14.1111 , we reproduce the linear- in- 
T correction (|2~52|) . 



5(T XX (T) ~ e 2 vV Tr 



(4.13) 



The constant term in (|4.11|) comes from the lower limit 
of the time integral, which is of the order of Ep . This 
constant merely renormalizes the bare value of the Drude 
conductivity. 

On the diffusive time scale (D)(D) ~ (£>£>), so that 
there is no difference between white-noise and smooth 
disorder. Therefore, in the diffusive limit the result (|4.1() 
applies, yielding the usual logarithmic correction 
In fact the contribution of the type 1)4. 1|) arises also in 
the ballistic regime when all terms in 12.38fl are taken 
into account. According to l|4.12|l . it has the form 



a xx D 



[ln(Tr) - const], 



(4.14) 



which is a subleading correction to the linear-in-T term 

(J233, EH- 

In the ballistic regime, Tr ^ 1, the above quali- 
tative arguments for a white-noise disorder can be re- 
formulated in terms of the interaction-induced renormal- 
ization of the differential scattering cross-section on a 
single impurity. Specifically, the renormalization occurs 
due to the interference of two waves, one scattered off 
the impurity and another scattered off the Friedel os- 
cillations created by the impuritjii&H. The interference 
contribution is proportional to the probability W(tt) of 
backscattering off the impurity (see Appendix and 
hence, to the return probability after a single-scattering 
event, as discussed above. 



On the other hand, this implies that the scatter- 
ing cross-section around <p ~ n is itself modified by 
the Friedel oscillations (in other words, the impurities 
are seen by electrons as composite scatterers with an 
anisotropic cross-section). The renormalization of the 
bare impurity depends on the energy of the scattered 
waves, which after the thermal averaging translates into 
the T-dependence of the effective transport scattering 
timeji& t(T) [this corresponds to setting t ~ T -1 in the 
return probability, see Eq. I|4.11|l ] . This mechanism pro- 
vides a systematic microscopic justification of the concept 
of temperature-dependent screening 17 . 

We recall that, in addition to the linear-in-T term, the 
conductivity correction contains a T-independent con- 
tribution determined by the ultraviolet frequency cut- 
off ~ Ep- In the case of strong interaction this term 
can be of the same order as the bare (non-interacting) 
Drude conductivity. The coefficient in front of this 
term cannot be calculated within the quasiclassical ap- 
proach because it is governed by short-distance physics 
at scales of the order of Xf- At the same time, ac- 
cording to the above picture, this T-independent correc- 
tion also modifies the impurity scattering cross-section 
around <j) = n. The corresponding correction dW(cj)) 
may thus be comparable to the bare isotropic scatter- 
ing probability Wq. An interesting consequence of this 
fact is a possible situation when the total relaxation rate 
t.T 1 cx J d0[Wo + (51^(0)] is smaller than the transport 
relaxation rate t~ x cx J d(f>[Wo + SW((f>)] (1 — cos<"6). 

In smooth disorder (small-angle scattering), the 
backscattering amplitude vanishes exponentially with 
kpd, and so does the amplitude of Friedel oscillations. 
This leads to the suppression of the Tr-contribution 
to the conductivity [see Sec. Ill CB : this fact was real- 
ized within the T-dependent screening picture already in 
Ref. [IT] for the case of scattering on long-range interface 
roughness]. We note, however, that the understanding of 
the interaction effects in terms of scattering off Friedel os- 
cillations is only possible in the ballistic regime. Indeed, 
the diffusive correction in a smooth random potential is 
not exponentially small and is related to random (having 
no 2fci?-oscillating structure) fluctuations of the electron 
density, as was pointed out in Refs. 0E1 The correla- 
tions of these fluctuations (which reduce to the Friedel 
oscillations on the ballistic scales) are described by the 
return probability at arbitrary scales. 

Finally, we use the interpretation of the interaction cor- 
rection in terms of return probability to estimate the MR 
in the white-noise random potential and at sufficiently 
weak magnetic fields, uj c <C T. Note that the zero-B bal- 
listic correction (|4.13|l does not imply any dependence 
of resistivity on magnetic field. Indeed, as follows from 
(|1.1|) . a temperature dependence of the transport time 
r(T) is not sufficient to induce any non-trivial MR, 

A Pxx (B,T) = p xx (B,T) - p xx (0,T) = 0, 

if r is ^-independent. 
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In order to obtain the independence of the resistiv- 
ity, we thus have to consider the influence of the mag- 
netic field on the return probability determining the cor- 
rection to the transport time. Since in the ballistic 
regime the characteristic length of relevant trajectories 
is L ~ vf/T -C I, their bending by the magnetic field 
modifies only slightly the return probability for lo c <C T. 
The relative correction to the return probability is thus 
of the order of (L/R c ) 2 ~ w 2 /T 2 independently of the 
relation between w c and r _1 . Therefore, to estimate the 
MR in the white-noise potential for oj c1 t _1 -C T, one can 
simply multiply the result (|4.13|) for B = by a factor 
{lu c /T) 2 , yielding 



Po 



Kt) 2 1 
kwl Tt' 



uj c < T. 



(4.15) 



A formal derivation of this result is presented in Sec. lVBl 
In a stronger magnetic field, lo c 3> T, the situation 
changes dramatically due to multiple cyclotron returns, 
see above. This regime is considered in Sec. IV Al below. 



MIXED DISORDER MODEL 



A. Strong B 



In Sec. IIIII we studied the interaction correction for a 
system with a small- angle scattering induced by smooth 
disorder with correlation length d 3> kp 1 . This is a typi- 
cal situation for high-mobility GaAs structures with suf- 
ficiently large spacer d. It is known, however, that with 
further increasing width of the spacer the large-angle 
scattering on residual impurities and interface roughness 
becomes important and limits the mobility. Furthermore, 
in Si-based structures the transport relaxation rate is 
usually governed by scattering on short-range impurities. 

This motivates us to analyze the situation when resis- 
tivity is predominantly due to large-angle scattering. We 
thus consider the following two-component model of dis- 
order ("mixed disorder"): white- noise random potential 
with a mean free time r wn and a smooth random po- 
tential with a transport relaxation time r sm and a single 
particle relaxation time r sniiS [r sm /r sm!S ~ (kpd) 2 1]. 
We will further assume that while the transport relax- 



ation rate r 1 



T^ 1 is governed by short-range 



disorder, r wn <C r sm , the damping of SdHO is dominated 
by smooth random potential, r smjS -C T wn . This allows 
us to consider the range of classically strong magnetic 
fields, uj c r wn 3> 1, neglecting at the same time Landau 
quantization (which is justified provided w c t SBIjS /7T <C 1). 

To calculate the interaction corrections, we have to find 
the corresponding kernel B a p(w,c\) determined by the 
classical dynamics. Naively, one could think that under 
the assumed condition T wn <C T sm the smooth disorder 
can simply be neglected in the expression for the classi- 
cal propagator. While this is true in diffusive limit, the 
situation in the ballistic regime is much more nontrivial. 



To demonstrate the problem, let us consider the kernel 
B X P J in the ballistic limit Tt ~ Tr wn 3> 1 and in a strong 
magnetic field lo c 3> T ^> t . If the smooth random po- 
tential is completely neglected in classical propagators, 
we have [see Appendix [BJ the second term in Eq. I|B.39|) 
can be neglected for uj c ^> T] 



B {p) 



1 



•%> , 2 1 

l frJ +9 ° 4\dQ 



( d Jl 



(5.1) 



where go(uJ, q) is the angle-averaged propagator with only 
out-scattering processes included, 



5o(w,q) 



in J^qRjJ-^qRc) 
u> c sin7r/i 



(5.2) 



and n = (uj + ir^)/uj c . If characteristic frequencies sat- 
isfy uj uj c (which is the case for T -C u> c ), Eq. I|5.2|) can 
be further simplified, 



9o 



4{Q) 

iuj + Twn 



(5.3) 



Substituting (|5.1|l and (|5.3|l in (|2.39|) . we see that 
momentum- and frequency-integrations decouple and 
that the first term in (|5.1(l generates a strongly 
ultraviolet-divergent g-integral ~ J dQ. 

The physical meaning of this divergency is quite trans- 
parent. The contribution of the first term in l|5.1|l to 5p xx 
is proportional to the time-integrated return probability 
J dtg (r = 0,t), similarly to 1)4. 1|) . For t -C r wn the 
propagator go (r, t) describes the ballistic motion in the 
absence of scattering, which is merely the undisturbed 
cyclotron rotation in the case of a strong magnetic field. 
Since at t — 2irn/uj c the particle returns exactly to the 
original point, the integral J dtgo(r = 0,t) diverges. 

The encountered divergency signals that the neglect of 
smooth disorder is not justified, even though r wn <C r sm . 
Indeed, with smooth disorder taken into account, the par- 
ticle does not return exactly to the original point after a 
cyclotron revolution, see Sec. IIVI The return probability 
is then described by Eqs. (|4.4|) . (|4.5|l with r replaced by 
r sm . It is worth mentioning a similarity with the prob- 
lem of memory effects in a system with strong scatterers, 
where even a weak smooth disorder turns out to be cru- 
cially importan1j2i£. 

To demonstrate the role of the smooth disorder on 
a more formal level, we write down the angle-averaged 
propagator in the ballistic regime, Tr wn ^> 1, for the 
mixed-disorder model, 



(2?( W ,q)) = 



4{Q) 



Q 2 /2r sm - iuj + t w 



(5.4) 



Clearly, in both limits r sm — oo and r wn = oo this for- 
mula reduces to l|5.3|) and (|3.2|l . respectively. In view of 
uJT wn 3> 1 the last term in the denominator of 1)5.4(1 can 
be neglected, and we return to the expression for solely 
smooth disorder. The presence of the term Q 2 /2r sm 
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regularizes the Q-integrals, thus solving the problem of 
ultraviolet-divergencies discussed above. The charac- 
teristic momenta are thus determined by Q 2 ~ Tr sm . 
Therefore, despite the weakness of the smooth disorder, 
r sm 3> r wn , it is the first (©/-dependent) rather than the 
third term which has to be retained in the denominator 
of 1)5.4(1 . In other words, in the ballistic regime and in 
a strong magnetic field the dynamics in the considered 
model is governed by smooth disorder. 

The above discussion demonstrates that at ui c T ^> 
t"^ the kernel B a p(u>, q) for the mixed-disorder model is 
given by (|2.38j) with propagators T> calculated in smooth 
random potential (i.e. with white-noise disorder ne- 
glected). The time r wn enters the result only through the 
matrices T a p (determined by the transport time r ~ r wn ) 
and S a p. Using r sm /r wn 3> 1, we find then that the re- 
sulting expression, 



B, 



1 



1 



2w 2 t w „ 



(VT>) 

[(V)(V)-2(Vn y )(n y V)} 



{(V)~2(n y Vn y )] 

c 

2 

— {nyDn x D) — (T>n x T>n x 'D) 



(5.5) 



is dominated by the first term corresponding to the first 
term in Eq. (|5.1|l . This yields for Q = qR c > 1 



W£T UO Z T (Q z — lU) Z 



4iQ) 



where f2 — 2ujT sm . 

As in previous sections, we first calculate the conduc- 
tivity correction for a point-like interaction. Substituting 
(G>3 in (232JI, we get 



Sa„ = -^uV (^) 



1/2 4cn 



(Tt)V2 : 



T > 1/t, 



(5.7) 

with the constant Co as defined in Eq. (|3.16(l . For an 
arbitrary (not necessarily small) value of the ratio t/t su1 
the coefficient 4 in (15.711 is replaced by 4— 3r/r sm . For r = 
r sm (i.e. without white-noise disorder) we then recover 
the ballistic asymptotics of Eq. (|3.15|l . 

As in the case of purely smooth disorder, the resis- 
tivity correction 5p xx is related to Sa xx via Eq. I|3.14[> . 
Comparing (|5.7|) with l|3.15fl . we see that the correction 
8p xx is enhanced in the mixed-disorder model by a factor 
~ 4(r sm /r) 1 ' /2 > 1 as compared to the purely smooth- 
disorder case. On the other hand, the scaling with tem- 
perature and magnetic field, Sp xx oc _B 2 T -1 / 2 , remains 
the same. 

Let us analyze now the crossover from the ballistic 
to the diffusive regime. Setting Tt ~ 1 in l|5.7|l . we 
find that the correction is parametrically large, 5a xx ~ 
(tsih/twii) 1 ^ 2 - Clearly, this does not match the diffusive 



contribution l|2.48[l . yielding 8a xx ~ 1 at Tt ~ 1, This 
indicates that returns without scattering on white-noise 
disorder continue to govern the correction in certain tem- 
perature window below T ~ 1/r, which normally belongs 
to the diffusive regime. 

To find the corresponding contribution, one should 
take into account the scattering-out term t~^ in the de- 
nominator of (|5.4|) . which yields 



1/2 



1 dzz^l 2 exp[— z/ttx] 



(2 7 ) 1/2 , 
4co7 1 / 2 x- 1 / 2 , 



sinh 2 ; 
x -C 1, 
X > 1, 



(5.8) 



where 7 = r sm /r 3> 1 and a; = Tt . To describe the 
temperature dependence of the interaction correction for 
all T, we have to add here the diffusive contribution, 
which has the form (|2.48|> for Tt <C 1 and vanishes for 
Tt ^> 1. This contribution corresponds to long times 
t ^S> t and describes the trajectories multiply scattered 
off white-noise disorder. Since at Tt ~ 1 the sum of the 
ballistic and diffusive contributions will be dominated by 
Gi(l,7) ~7 1/2 > 1, the precise way of vanishing of the 
diffusive contribution at Tt ~ 1 is inessential. Therefore, 
we can describe it by the function Gq(x), Eq. I|3.12[) . The 
resistivity correction for a system with mixed disorder 
and point-like interaction has thus the following form: 



5pxx{B) (w c t) s 



Po 



vV Gf x (TT,T sm /T), (5.9) 



(5.6) where 



GT x (x, 7 ) 



G 1 (x,j) + G (x) 

' -lnx + (27) 1 / 2 , 
4c o7 1 / 2 x- 1 / 2 , 



(5.10) 



x < 1, 
x > 1, 



This result is illustrated in Fig. II lb . 

In the case of Coulomb interaction, we have as usual a 
similar result for the exchange contribution 



Spl^(B) ( Wc r) ; 



po 



ixkpl 



Gf X (TT,T sm /T), 



(5.11) 



with 



GJ? IX (x, 7 ) 



-Gi(x, 7 ) + G F (x) 
f -lna;+( 7 /2) 1 / 2 , 



(5.12) 



2c 7 



1/2-1/2 



x < 1, 
x> 1. 



This function is shown in in Fig. II lb . In fact, here the 
diffusive contribution can be described either by the func- 
tion Gf (x) or by Go (x) because in the diffusive limit they 
coincide up to a small constant. Since in the intermedi- 
ate and ballistic regimes [where Gp(x) and Gq(x) differ] 
the contribution ^G%(x, 7) is dominant, the behavior of 
the diffusive contribution is of no importance, as in the 
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\Fq /(l + Fq ) in the ballistic term. The result thus reads 



FIG. 11: Functions G™ ix (Tr) (a) and G™ ix (Tr) (b) describ- 
ing the temperature dependence of the resistivity correction 
due to point-like and Coulomb (exchange) interaction, respec- 
tively, in the mixed-disorder model for different values of pa- 
rameter 7 =. r sm /r = 20, 10, 5 (from top to bottom). Dashed 
curves represent these functions for purely smooth disorder 
( 7 = 1). 



case of point-like interaction. Note that the ballistic con- 
tribution corresponds to the point-like interaction with 
vVq = i, yielding a factor i in front of Gi(x,~f) as com- 
pared to H5.10fl . This is because the dynamical part of 
screening is suppressed for all relevant Q ~ Tr sm 3> 1 in 
the whole range of temperatures, even for Tr < 1, where 
this contribution is important. 

This also applies to the Hartree contribution to the 
resistivity. Within the u Fq- approximation" we have 
again an effectively point-like interaction with vVq ~ 



^ mix (£) 

PO 

where 

Gfr(*,7) = 
1 



= 3^-^G£ IX (Tr,r sm /r), 



(5.13) 



F 



= < 



2 (1 + Fg) 
ln(l + Fg) 



°— G 1 (x, 7) + G\ i (x) 



-2cn 



(1 + F C 



In a; 



-7 1 / 2 *- 1 / 2 , 



7 \ 1/2 

2, 



x < 1, 

X > 1. 

(5.14) 



Before closing this subsection, we briefly discuss the 
Hall resistivity in the mixed disorder model. Repeating 
the steps described above, we find that the ballistic con- 
tribution to p xy also contains an extra factor (rgm/r) 1 / 2 , 
similarly to p xx . For an arbitrary (not necessarily small) 
value of the ratio T/r sm the coefficient 11 in Eqs. I|3.52[l . 
(pTBljl is replaced by [6 + 5r/r sm ] (rsm/r) 1 / 2 . 



B. Weak B 

In the case of a purely smooth disorder (Sec. IIII|) the 
resistivity correction in the ballistic regime is exponen- 
tially suppressed for oj c <^ T because the particle cannot 
return to the origin. When the short-range potential is 
present, the situation changes and the return probability 
is determined for Tr ^> 1 by the single-backscattering 
processes. The interaction-induced MR arises then due 
to the influence of the magnetic field on the probabil- 
ity of such return, as discussed in the end of Sec. IIVI 
In this case, there is no need to take the smooth poten- 
tial into account and the MR is determined solely by the 
white-noise disorder. Let us calculate the corresponding 
correction using the ballistic form JB 



of the kernel 



ABjx ■ 

For a point-like interaction, substituting l|B.40|) in 
(|2.39|l . we find the following ballistic (Tr ^> 1) asymp- 
totics of the longitudinal MR, 



po 



(uj c t) 2 vVq it 
irkpl 727V' 



(5.15) 



In the case of Coulomb interaction, AB X P J is multiplied 
by the ballistic asymptotics of the interaction, Eq. (|C.ip . 
Substituting this product in Eq. (|2.39l) . we get the Fock 
contribution to the MR in the form 



XX 



Po 



(lo c t) 2 17tt 
irk F l 192 Tr' 



Tr > 1. 



(5.16) 



The corresponding Hartree term also scales as B 2 /T. It 
is worth noting that there is another contribution to the 
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MR in this regime, which comes from the suppression 
of the triplet channel due to Zeeman splitting E% rather 
than from the orbital effects. This contribution is iden- 
tical to that found in Ref. J2(J for the ballistic magne- 
toresistance in a parallel magnetic field. It also scales 
as B 2 JT in a weak magnetic field; however, it contains 
an extra factor (Ez/uj c ) 2 , as compared to <|5.16|l . This 
factor is small in typical experiments on semiconductor 
heterostructures where the effective mass of the carriers 
is much smaller that the bare electron mass. 

We now turn to the Hall resistivity. Using IB. 381) and 
(|B.41|I . we find for oj c , r _1 <C T and for arbitrary relation 
between to r and r _1 



Pxy 



Trkpl 12 TV 



for the point-like interaction, and 
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xy 



Pxy 



1 - 



49(w c r) 2 



330 



llTT 

96 Tt 



(5.17) 



(5.18) 



for the Coulomb interaction. The result 1|5.18|) reduces 
in the limit B — ► to that obtained in Ref. [RJ from the 
quantum kinetic equation. We see that in view of a rela- 
tively small value of the numerical coefficient 49/330, the 
first (B-independent) term in square brackets in (I5.18|) 
dominates for uj c t < 1, so that the results of Ref. [Rj are 
applicable in sufficiently broad range of magnetic fields. 
For the corresponding Hartree-correction to 8p xy calcu- 
lated within the "i^-approximation" , we refer the reader 
to Ref. E 



VI. ANISOTROPIC SYSTEMS 
A. Qualitative discussion 

In the preceding consideration, we assumed that the 
2D system is isotropic. While this is true for the majority 
of magnetotransport experiments we have in mind, there 
exists a number of important situations when the trans- 
port is anisotropic, a xx ^ a yy . First, such an anisotropy 
can be induced by the orientation of the 2D electron gas 
plane with respect to the crystal axes, see e.g. Ref. |48| 
for a measurement of the quantum correction for the 
(110) surface of the Si-MOSFET. Second, the electron- 
electron interaction may lead to spontaneous formation 
of a charge-density wave. Finally, the anisotropy may be 
induced by a one-dimensional periodic modulation (lat- 
eral superlattice). The latter example is of special inter- 
est in view of emergence of commensurability oscillations 
(known as Weiss oscillations)^, and we will discuss it in 
more detail in Sec. IVI CI 

The interaction-induced correction to the conductiv- 
ity tensor of an anisotropic system was calculated for 
the diffusive regime and B = by Bhatt, Woffle, and 
Ramakrishnanii. They showed, in particular, that the 



quantum correction preserves the anisotropy of the qua- 
siclassical (Boltzmann) conductivity. Below we will gen- 
eralize their result onto the case of a classically strong 
magnetic field, and, furthermore, will extend the consid- 
eration to the ballistic regime. 

We begin by presenting a simple argument allowing one 
to estimate the conductivity correction in an anisotropic 
system; we will confirm it by a formal calculation in 
Sec. IVI Bl According to Eq. (|4.1(l . the relative correction 
to a diagonal component er^ (/i = x, y) of the conductiv- 
ity tensor is determined by the return probability (and 
is, thus, the same for [i — x and ji = y). This implies, in 
the diffusive regime 



5a u 



-Re 



1 



yielding 



(dq) 



1 



B a jiq a qp 



u=l/r 



uj=T 



Sa x 



1/2 



In Tt 



(6.1) 



(6.2) 



and analogously for Sa yy . In the ballistic regime the 

rp-l 

time-integrated return probability J dt(D(t)) scales 
as (Tt)~ 1/2 [see Eqs. (|4.5jl and l|4.6(l ]. so that we have 
instead of l|6.2|l . 



6a xx ~ e z K — (Tt) 



\-l/2 



(6.3) 



The explicit form of the function JC(x) will be calculated 
below. Since the conductivity corrections (|6.2|) and l|6.3|) 
are only determined by the anisotropic diffusion, we ex- 
pect that they do not depend on the particular source 
of anisotropy, in analogy with Ref. l4ll An important 
feature of the results 1)6. 2|) and (|6.3|l is that they mix the 
components a xx and a yy of the conductivity tensor. This 
will play a central role in our analysis of the interaction 
effect on the magnetoresistivity of modulated systems in 
Sec. lvTTl 

It is worth mentioning that the validity of the formula 
(|6.3() for the ballistic regime is restricted on the high- 
temperature side by the condition T < T a d, where T~^ 
is the time scale on which the anisotropic diffusion of the 
guiding center sets in. The value of T a( j depends on the 
particular microscopic mechanism of anisotropy. We will 
estimate T a( j and the behavior of the conductivity cor- 
rection at T ^> T a( j for a modulated system in Sec. IVI 01 



B. Calculation of the interaction-induced 
correction to resistivity 

We proceed now with a formal calculation of the quan- 
tum correction to the conductivity of an anisotropic sys- 
tem in a strong magnetic field. As a model of anisotropy, 
we will assume anisotropic impurity scattering, with a 
cross-section W(4>,4>') ^ W(4> — (/)'). Repeating the 
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derivation performed in Sees. Ill Al and ITTbI we find that 
the result 1)2. 35[) . (|2.36l) remains valid in the anisotropic 
case, with the matrix T a p proportional to the correspond- 
ing (anisotropic) diffusion tensor D a 0, 



T. 



a/3 



2D a p 



1 + UJ?T x T y 



^c^~x^~y 



-U) c T x T y 
T,, 



(6.4) 

where t x and t v are the relaxation times for the corre- 
sponding components of the momentum. We begin by 
considering the diffusive limit, when the leading contri- 
bution comes from three-diffuson diagrams, Fig.^? and e 
(sec Sec. Ill Cl O. which are represented by the last term in 
Eq. I|2.36[) . The singular contribution to the propagator 
T>, governed by the diffusion mode, has a form analogous 

to ifmji . 

V (w, q; (p, <P ) ^ =: ; , (6.5) 

D a pq a qf3 - iu 

see Appendix [E] f° r the derivation of l|6.5|) and explicit 
expressions of <5 r,l- Using l|6.5[) and l|E.3|) . we get 



(V) ~ (V s ) 



1 



D al 3q a qf3 - ico 



(6.6) 



and 



B M (w,q) ~ ~(V s n a V s n p V s ) 
4 D 2 q 2 

XX^X 



v 2 F (D xx q 2 + D yy q 2 - iuj) 3 



(6.7) 



The result (|6.7|> can also be obtained with making use 
of the identity J33j; then it is sufficient to keep only 
the leading term (unity) in the expressions for functions 
^r ,l en tering ijfTH)) . Substituting l|67r)) . gSSJl, lEH 
in H2.35|l . we obtain the final result for the conductivity 
correction in the diffusive regime, 



Sa x 



2tt 2 \D 



D, 



1/2 



In Tr, 



2tt 2 \D 



D 



1/2 



In Tr, 



(6.8) 
(6.9) 



in full agreement with a qualitative consideration of 
Sec. rvTAl [Eq. (|6~2"|) ]. The correction to the Hall con- 
ductivity is zero in the leading (In Tr) order, as in the 
isotropic case. For the point-like interaction, the result 
remains the same, up to a factor vVq. 

We now extend the consideration beyond the diffu- 
sive limit (thus allowing for qR c > 1), assuming first the 
smooth disorder and concentrating on longitudinal com- 
ponents of the conductivity and resistivity tensors. In 
analogy with JOJ, the singular contribution V s to the 
propagator acquires then the form (see Appendix lE|) 

V(lo, q; (/), (j>') = exp{-iqR c [sin(4> - cj> q ) - sin(0' - (f> q )]} 

„ x{<t>)xi<t>') 



D a/3 q a qp - iu) 



(6.10) 



where 

iqvF ( 1 1 

X\4>) = 1 T~ cos 4> cos <t>q H sm H> sm < 



This yields 



(V) 



J$(qR c ) 



and 



B xx (u>,q) = 



D xx ql + D yy q 2 - iu 

4 J 2 (gR c )D 2 xx q 2 x 
v% (D xx q 2 + D yy q 2 y - iu;) 3 ' 



(6.11) 
(6.12) 

(6.13) 



which differs from l|6.6|) . I|6. 7fl by the factor Jq (qR c ) only. 
In the ballistic limit Tt x ,Tt v ^> 1 the relevant values of 
qR c are large, qR c 3> 1, so that the screening is effec- 
tively static and the interaction is effectively point-like 
with Vq = l/2v. Substituting then ljfi.13^ in (|2.35|) and 
rescaling the integration variables q x = D xx ^ 2 q x , q y = 

— 1/2 

D yy q y , we find 



8a 



ii y 



^o(Tr y )^ 2 



x -K \/ 1 — D xx / D 



7T 

c ^yy c 

0CT xx — "TT OCT yy , 

'J XT 

where K is the elliptic integral, 
dx 



yy I > 



(6.14) 
(6.15) 



cos z x + q z sm x 



K(Vl-9 2 ), 0<?<1, 
(6.16) 

and we assumed that y is the easy-diffusion axis, D yy > 

D X x • 

Let us analyze the obtained results in the limits of weak 
and strong anisotropy. It is convenient to set D xx = -Do, 
t x = To, D yy = Dq + AD, and to introduce a dimen- 
sionless anisotropy parameter a = AD/Dq. Using the 
asymptotics of the elliptic integral, 



K(s) 




we find 



s < 1, 
1 - s « 1, 



8a xx ^ -^2 c o( TT y) 1/2 

l--r, a<l, 



(6.17) 



ln(16a) 

7TO? 1 / 2 



(6.18) 



a > 1, 



and 8a yy = (1 + a)8a xx Equations (|6.14() . H6.15|I and 
(|6.18|l confirm the qualitative arguments in Sec. IVI Al 
(based on the behavior of the return probability) which 
led to Eq. iftTSjl . 



24 



C. Modulated systems 



In this subsection, we apply the results of Sec. lVIBl to 
a particularly important class of anisotropic 2D systems, 
namely, 2D electron gas subject to a periodic potential 
varying in one direction. Such systems (lateral superlat- 
tices) have been intensively investigated experimentally 
during the last fifteen years. In a pioneering work 49 , 
Weiss, von Klitzing, Ploog and Weimann discovered that 
even a weak one-dimensional periodic modulation with 
a wave vector k || may induce strong oscillations of 
the magnetorcsistivity p xx (B) [while showing almost no 
effect on p yy (B) and p xy (B)) 1 with the minima satisfy- 
ing the condition 2R c /a = n — 1/4. Here n = 1,2,... 
and a = 2n/k is the modulation period. The quasiclas- 
sical nature of these commensurability oscillations was 
demonstrated by Beenakker— , who showed that the in- 
terplay of the cyclotron motion and the superlattice po- 
tential induces a drift of the guiding center along y axis, 
with an amplitude squared oscillating as cos 2 (kR c — 7r/4) 
[this is also re produced by a quantum-mechanical calcula- 
tion, see Ref s. 51 52 53]. While Ref . l5(J assumed isotropic 
impurity scattering (white-noise disorder), it was shown 
later that the character of impurity scattering affects cru- 
cially the dependence of the oscillation amplitude on the 
magnetic field. The theory of commensurability oscilla- 
tions in the situation of smooth disorder characteristic for 
high-mobility 2D electron gas was worked out in Ref. |w| 
(see also numerical solution of the Boltzmann equation 
in Ref. I55T ) and provided a quantitative description of 
the experimentally observed oscillatory magnetoresistiv- 
ity Ap xx (B). The result has the form^i 



Ap xx WT] k IR, 



po 4 sinh(7rA) 



J tX {kR c )J- lX (kR c ), (6.19) 



where r/ is the dimensionless amplitude of the modulation 
potential \V(x) — r\Ep cos(fcx)], and 



A = 



1 



1 - 



\ + -{kR c y 

T 



"1/2 



(6.20) 



In the range of sufficiently strong magnetic fields 
Eq. I|6.19|) describes the commensurability oscillations 
with an amplitude proportional to B 3 , 



Ap x 



Po 



irkRc 



cos 2 (kR c -7r/4), 



(6.21) 



For precise conditions of validity of Ij6.21|l , as well as for 
an analysis of the result l|6.19fl in the whole range of 
magnetic fields, the reader is referred to Ref. 54. 

As to the modulation-induced corrections Ap xy , Ap yy 
to the other components of the resistivity tensor, they are 
exactly zero within the quasiclassical (Boltzmann equa- 
tion) approach, independently of the form of the impurity 
collision integralS^i^. Such corrections appear in the 
quantum-mechanical treatment of the problem£§i££ and 
are related to the de Haas- van Alphen oscillations of the 



density of states induced by the Landau quantization of 
spectrum. As a consequence, these oscillations are expo- 
nentially damped by disorder, with the damping factor 
~ exp[— 2it/uj c t s ]. The phase of such quantum oscilla- 
tions A/9yy OS ' ) is opposite to that of quasiclassical com- 
mensurability oscillations in Ap xx , Eqs. (|6.19() . 16.21f) . 
Indeed, oscillations in Ap yy that are much weaker than 
those in Ap xx , have the opposite phase, and vanish much 
faster with decreasing £?, were observed in Ref. We 
will neglect these oscillations, which are exponentially 
weak in the range of magnetic fields considered in the 
present paper, uj c t s /it <C 1. We are going to show that 
the interaction-induced correction to resistivity also gen- 
erates oscillations in p yyj which are, however, unrelated 
to the DOS oscillations of a non-interacting system and 
become dominant with lowering temperature. 

To demonstrate this, we apply the result of Sec. IVI Bl 
for the interaction-induced correction in an anisotropic 
system. The anisotropy parameter is governed by the 
quasiclassical correction to p xx due to modulation, 



® xx 



1 



p xx 
Pyy 



1 = 



Ap xx 
Po 



(6.22) 



and is given by Eq. I|6.19fl . For simplicity we will as- 
sume that the effect of modulation is relatively weak, 
a -C 1. (Generalization to the large-a case with making 
use of the corresponding formulas of Sec. IVf Bl is com- 
pletely straightforward.) Using (|6.8|) and l|6.15|l . we find 
the oscillatory correction to p yy as a combined effect of 
the modulation and the Coulomb interaction, 



Sp 



in) 



po 



(uj c t) 2 Ap xx 
lixkpl po 
- In Tr, 



1/2 



Tt< 1, 
Tr> 1. 



(6.23) 



In the presence of strong scatterers (mixed disorder 
model), the result for the ballistic regime is enhanced 
by the factor 4(r sm /r) 1 / 2 ^> 1, as discussed in Sec.lVl 

Let us remind the reader that the result i|ti.23f) is valid 
for temperatures T <C T a d, where T", 1 <C r is the char- 
acteristic time on which the motion of the guiding center 
takes the form of anisotropic diffusion (see Sec. IVI A|) . 
For the case of a modulated system with a smooth ran- 
dom potential we find T~~r ~ r(a/i? c ) 2 . This is because 
on a scale shorter than the modulation leads to a 
drift of the guiding center along y axis with the velocity 
depending on the coordinate X of the guiding center—, 

v d (X) = --^kR c J (kR c )sm(kX) (6.24) 



rjVF 



y/2^kR t 



cos(kR c - 7r/4) sm(kX). 



In a time a 2 /D = the position X of the guiding- 
center is shifted by a distance of the order of the mod- 
ulation period a due to the small-angle impurity scat- 
tering. Therefore, the drift velocity Vd typically changes 
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sign on this time scale, so that the drift is transformed to 
an additional diffusion process, with AD yy ~ (^)T, 



ad ' 



in agreement with (j6.21|) . To estimate the resistivity cor- 
rection Spyy in the ultra-ballistic regime T 3> T a d, we use 
the relation between the conductivity correction and the 
return probability fSec. IfVjl . The return probability R„ 
after n revolutions (introduced in Sec. HV|l is modified by 
the modulation-induced drift in the following way: 



R n 1 - n U C T 



2v\ 



(6.25) 



According to (|4.1|) , this yields an oscillatory correction to 
resistivity suppressed by a factor ~ T a d/T as compared 
to the second line (ballistic regime) of Eq. H6.23|l . 

Let us summarize the results obtained in this sub- 
section. We have shown that in a periodically modu- 
lated system the interaction induces, in addition to the 
quadratic MR studied in Sees. IIIII and IIVI an oscillatory 
contribution to the component p yy of the resistivity ten- 
sor, which is not affected by modulation (and thus shows 
no oscillations) within the Boltzmann theory. When the 
parabolic MR is negative (meaning that the exchange 
contribution dominates), which is the case under typical 
experimental conditions and for not too high tempera- 
tures, these quantum interaction-induced oscillations in 
Pyy are in phase with classical oscillations in p xx , as 
follows immediately from Eq. I|6.23f) [see Fig. In 
other words, their phase is opposite to that of the above- 
mentioned contribution Apy^ OS ^ induced by the DOS os- 
cillations. 

We come therefore to the following conclusion con- 
cerning the phase of the total oscillatory contribution 
to pyy. While at sufficiently high temperatures the p yy - 



(DOS) 



[and 



oscillations have, due to the contribution Ap 
possible due to the Hartree counterpart of Eq. 
the phase opposite to Ap xx , with lowering temperature 
the exchange contribution Eq. i|tj.23f) starts to dominate, 
implying that p yy oscillates in phase with p xx . Fur- 
thermore, the both contributions are damped differently 
by disorder: the high-temperature out-of-phase oscilla- 



tions Ap 



(DOS) 



vanish with lowering B much faster that 



the low-temperature in-phase interaction-induced oscil- 
lations 5 Pyy. 

Our results are in qualitative agreement with a recent 
experiment^. It was observed there that at sufficiently 
high temperature, T > 2.5 K 1 the oscillations in p yy 
have the opposite phase with respect to p xx , in accord 
with earlier experimental findings 49 . However, when the 
temperature was lowered, the phase has changed and p yy 
started to oscillate in phase with p xx , with an amplitude 
increasing with decreasing T. In addition to these novel 
oscillations, a smooth negative MR was seen to develop 
in the same temperature range. The authors of Ref. |58| 
emphasized a puzzling character of the temperature de- 
pendence of the observed oscillations, which cannot be 
explained by earlier theories£2i£ii£2i£2iii4 discarding the 
interaction effects. Our theory leading to Eq. l|tj.23fl pro- 
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FIG. 12: Magnetoresistivity in a lateral superlattice with 
modulation wave vector k || e^. (a) Quasiclassical Weiss- 
oscillations; dashed line shows the resistivity in the absence 
of modulation, (b) Interaction-induced quantum oscillations 
in p yy for three temperatures. The curves correspond to the 
values of the parameter 2co(r sm /Tr 2 ) 1 ^ 2 = 0.1, 0.3, 0.5 (from 
top to bottom), assuming mixed disorder. Dashed line rep- 
resents the resistivity of the non-interacting system. Typi- 
cal experimental parameters are used: effective mass m = 
0.067 x 9.1 x 10" 28 g, electron density n e = 3.16 x 10 11 cm -2 , 
modulation strength n = 0.05, modulation period a — 260nm, 
momentum and single-particle relaxation times r = lOOpc 
and t s = 5pc, respectively. 



vides a plausible explanation of these experimental find- 
ings. Quantitative comparison of the theory and experi- 
ment requires, however, a more systematic experimental 
study of the temperature dependence of the amplitude of 
Pyy-oscillations in a broader temperature range. 
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VII. CONCLUSIONS 
A. Summary of main results 

Let us summarize the key results of the present pa- 
per. We have derived a general formula Ij2.35|) . (|2.38() for 
the interaction-induced quantum correction 8a a p to the 
conductivity tensor of 2D electrons valid for arbitrary 
temperature, magnetic field and disorder range. It ex- 
presses 5a a {3 in terms of classical propagators in random 
potential ( "ballistic diffusons" ) . In the appropriate lim- 
iting cases, it reproduces all previously known results on 
the interaction correction (see Sec. Ill C|) . 

Applying this formalism, we have calculated the inter- 
action contribution to the MR in strong B in systems 
with various types of disorder and for arbitrary Tt. In 
the diffusive limit, Tt <C 1, the result does not depend on 
the type of disorder, as expected. Specifically, the MR 
scales with magnetic field and temperature as follows, 
Spxx oc B 2 ln(TV) and Sp xy oc B In Tt, in agreement with 
Refs.[l3and[lll 

In the ballistic limit, Tt ^> 1, the result is strongly 
affected by the character of disorder. In Sec. lIIII we have 
performed a detailed study of the case of smooth dis- 
order characteristic for high-mobility GaAs heterostruc- 
tures. We have found that the temperature-dependent 
MR scales at u c > T as Sp xx oc B 2 (Tt)~ 1/2 and 
5p xy oc B(Tr) 1 / 2 . In addition, there is a temperature- 
independent (but larger) contribution oc B 3 / 2 to the Hall 
resistivity. In the opposite limit lo c <C T the MR is sup- 
pressed. 

We have further considered a mixed disorder model, 
with strong scatterers (modeled by white-noise disorder) 
superimposed on a smooth random potential (Sec. EJ- A 
qualitatively new situation arises when the momentum 
relaxation rate r^ 1 due to smooth disorder is much less 
than the total momentum relaxation rate t _1 . Such a 
model is believed to be relevant to Si-based structures, 
as well as to GaAs structures with very large spacer. We 
have shown that in the ballistic limit and at lu c 3> T the 
corrections to both longitudinal and Hall resistivities are 
enhanced (as compared to the case of smooth disorder) 
by a factor ~ (Tsm/r) 1 / 2 ^> 1. In the range of weaker 
magnetic fields, oj c <C T, the interaction-induced MR 
scales in the ballistic regime as Ap xx oc i? 2 (Tr)~ 1 and 
5p xy oc _B(Tr) _1 [l — const (lu c t) 2 ]. 

For a weak interaction (k -C hp) the correction is dom- 
inated by the exchange contribution, implying that Ap xx 
is negative and that the slope of p xy decreases with in- 
creasing temperature. This is true up to a tempera- 
ture Th 3> t -1 (defined in Sec. IIII D|) where the sign 
changes. In the case of a strong interaction the mag- 
nitude of the Hartree contribution (and thus the sign 
of the total correction) depends on angular harmonics 
FP' a of the Fermi-liquid interaction (Sec. IIIIE|) . It is 
worth emphasizing that in contrast to the diffusive limit 
where only Fq is relevant, in the ballistic regime all the 
Fermi-liquid parameters are, strictly speaking, impor- 



tant, see Eq. i|3.40|l . Therefore, predictions of the u Fq- 
approximation" (with only one Fermi-liquid parameter 
retained) should be treated with caution. 

We have further applied our formalism to anisotropic 
systems (Sec. IVIfl and demonstrated that the correction 
mixes the components p xx and p yy of the resistivity ten- 
sor. This result is of special interest in the case of systems 
subject to a one-dimensional periodic modulation (lateral 
superlattice; wave vector k || e x ). Specifically, we have 
shown that the interaction induces novel oscillations in 
p yy , which are in phase with quasiclassical commensura- 
bility (Weiss) oscillations in p xx . 



B. Comparison with experiment 

Our results for p xx in the case of smooth disorder (pub- 
lished in a brief form in the Letter—) have been confirmed 
by a recent experiment on n-GaAs system^, which was 
performed in the broad temperature range, from the dif- 
fusive to the ballistic regime. Specifically, Li et a/— 
found that the MR scales as Ap xx oc B 2 in strong mag- 
netic fields. The obtained temperature-dependence of the 
proportionality coefficient G(Tt) was in good agreement 
with our predictions. 

Very recently, Olshanetsky et a/— studied the MR 
in the ballistic regime in a Si/SiGe structure of n-type, 
where both short- and long-range potential are expected 
to be present. They found the interaction-induced cor- 
rection to p xx larger by a factor ~ 5 as compared to 
our prediction^ for the case of smooth disorder. This 
conforms with the results of the present paper for the 
mixed-disorder model, where we find an enhancement of 
Ap xx by a factor 4(r sm /r) 1 / 2 > 1. 

As has been mentioned in Introduction, the 
interaction-induced MR in the ballistic regime was mea- 
sured for the first time as early as in 1983, by Paala- 
nen, Tsui, and HwangiS, who studied GaAs structures. 
Again, a parabolic temperature-dependent MR Ap xx was 
obtained, in agreement with our findings. However, its 
magnitude was considerably (roughly an order of mag- 
nitude) larger compared with our theoretical result for 
the case of smooth disorder, as well as with the re- 
cent experiment. We speculate that samples used in 
Ref. ^2 probably contained an appreciable concentration 
of background impurities, which has led to an enhance- 
ment of the interaction-induced contribution to resistiv- 
ity, similarly to the recent work— (Indeed, the results 
for the mixed-disorder model shown in Fig. 1111 may cre- 
ate an impression that the log T behavior extends up to 
Tt ~ 10, as was concluded in Ref. [12J.) Remarkably, the 
interaction-induced quantum correction to conductivity 
may serve as an indicator of the dominant type of disor- 
der. 

To the best of our knowledge, no experimental study 
of the interaction effects on Hall resistivity p xy has been 
published. This part of our predictions therefore awaits 
its experimental verification. 



27 



Finally, our results for systems with one-dimensional 
periodic modulation are in qualitative agreement with 
the recent work by Mitzkus et al,^, as we discussed in 
detail in Sec. IVI CI Quantitative comparison of the the- 
ory and experiment requires an experimental study of the 
temperature-dependence of novel oscillations (found ex- 
perimentally in Ref. and theoretically in the present 
paper) in a broader temperature range. 



Outlook 




VvA/W 



FIG. 13: Aslamazov-Larkin-type diagrams describing the 
Coulomb-drag contribution to the resistivity, which cancels 
the "inelastic" part of the diagrams /, g from Fig. 



Before closing the paper, we list a few further ap- 
plications of our formalism and its possible generaliza- 
tions. First, our results can be generalized to frequency- 
dependent (rather than temperature-dependent) MR. 
Second, the interaction effects in systems of other di- 
mensionality, as well as in macroscopically inhomo- 
geneous systems, can be investigated by our general 
method. Third, the formalism can be used to calculate 
the phonon-induced contribution to resistivity, which be- 
comes larger than that due to Coulomb interaction at 
sufficiently high temperatures. Further, thermoelectric 
phenomena in the full range of magnetic fields and tem- 
peratures can be studied in a similar way. Finally, our 
approach can be generalized to the regime of still stronger 
magnetic fields, where the Landau quantization can not 
be neglected anymore; the work in this direction is in 



progres: 
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APPENDIX A: CANCELLATION OF THE 
INELASTIC TERM 

As discussed in Sec. Ill Al diagrams / and g give rise, in 
addition to the contribution (|2.25() , to a term of the type 
()2.26(l . characteristic for inelastic effects. This yields at 
B = a disorder-independent correction to resistivity 
5p ~ (T/eEp) 2 , see below. Note that such a contribu- 
tion to resistivity would be obtained if one substitutes 
the inelastic relaxation rate of a clean 2D electron gas, 
Tj^ el ~ T 2 /Ef in the Drude formula However, such 

a procedure clearly makes no sense. Indeed, in a transla- 
tionally invariant system electron-electron collisions con- 
serve the total momentum and thus give no contribution 
to resistivity. Therefore, the correction l|2.26|l ought to 
be canceled by some other contribution. Below we show 



explicitly that this is indeed the case, and that this sec- 
ond contribution is of the Coulomb-drag type, described 
by the diagrams in Fig. IT51 

For simplicity, we restrict our consideration here to the 
case of zero B and white-noise disorder, which allows us 
to use the results of Ref. |6l| for the Coulomb drag. Note 
that while Ref. considered the drag between two lay- 
ers, we refer to the "self-drag" within a single layer. As 
we will see below, the characteristic momenta q determin- 
ing the contribution (|2.26|) are large, q ~ kp. For this 
reason, there is no need to take into account impurity-line 
ladders while evaluating this term, similarly to the calcu- 
lation of drag in Ref. for a small inter-layer distance. 
We thus have 



1 



2'KVV 2 F 



d 2 p 

(2n) 2 



x Re [2p 2 x G 2 l {e 1 p)G K {e -ui,p- q)G A (e,p) 

+ Px(p x - q x )G R {e,p)G A {e,p) 

x G R (e-cj,p-q)G A (e-u,p-q)], (A.l) 

where G RjA (e,p) = (Ep + e — p 2 /2m ± i/2r) _1 are the 
disorder-averaged retarded and advanced Green's func- 
tions. Using the identity 

G R (e,p)G A (e,p) = ir[ G R (e,p) - G A {e,p) ], 

we reduce l|A.l(l to the form 

q) = -^/ tJ^p^ 

x Re[G R (e,p)G A (e-cj,p-q)] 
- -^-Iml%,q), (A.2) 



2v F i/cu 



where n(o;,q) is the polarization operator l|2.15[) . 
f d 2 

Imn(w,q) = - / -— 1 G K {t,p)G A {e-uj,p-q) 



it J (2tt) ; 



it) 



2v 6(qvF — uj) 

qvp 



(A.3) 



where 9(x) is the step function. Furthermore, the imag- 
inary part of the interaction propagator within the RPA 
is proportional to ImII(a;,q) 

lmU(uj,q) = -|[/(w,q)| 2 Imn(w,q). (A.4) 
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Substituting (|A.2|I and (|A.4|I in (|2.26(l , we finally obtain 
,oa d(J 1 f d 2 q 



5a™ 1 = 



27r2Tsinh 2 (w/2T) J (2tt) 



x g 2 |C/(w,q)| 2 [Imn(u;,q)] 2 , 



(A.5) 



This expression is identical, up to a sign, to the result of 
Ref. for Coulomb drag. This demonstrates that two 
contributions indeed cancel each other, 



5a" 



<5cr drag = o. 



(A.6) 



Using the explicit form of lmH(cj,q), Eq. (|A.3|I . and of 
U(oj,q), Eq. H2.14J) . in the ballistic regime, it is easy to 
estimate Sa lnel (we assume here k ~ &f for simplicity), 



5a" 



-e 2 T 2 T 



did 



^ ^ ~ -e 2 (Tr) 2 . (A.7) 



As has been stated above, the g-integral is determined 
by the ultraviolet cutoff. 

Finally, we note that in double-layer system the inter- 
layer interaction does give rise to a correction 5a lncl to the 
driving-layer conductivity, which is equal in magnitude 
and opposite in sign to the transconductivity. This effect 
is, however, reduced by a factor ~ (fc_F^)~ 4 (where £ is the 
interlayer distance), as compared to (|A.7() . see Ref. loll 



APPENDIX B: PROPAGATOR AND KERNELS 
B a/3 FOR WHITE-NOISE DISORDER 

In this appendix we will derive the general expressions 
(valid for arbitrary magnetic field) for the kernels B X P J 
and Bx P y in terms of the quasiclassical propagator for a 
white-noise random potential. This will allow us to re- 
produce the results of Refs. lTMl^ . where the interaction- 
induced corrections to a xx and p xy were studied for a 
white-noise disorder in the limit B — > 0. We will further 
apply the formalism to calculate the longitudinal MR 
and the Hall resistivity in a finite magnetic field with 
lo c <C T. The resistivity tensor in yet stronger magnetic 
field, oj c 3> T, is studied, in the more general framework 
of a mixed disorder model in Sec. [Vj 

Using Eqs. and we get 

B x p } = ±(V) 2 -±(Vn x )(n x V) 
It t 

1 



+ -(V)-(n x Vn x ) 



~(n x T>n x T>) + 2u} c (n x Vn y V) 
., (Vn x T>n x T>) 



1-u 2 t 2 



+ ^-{Vn x Vn y V) 



(B.l) 



for the kernel describing the longitudinal resistivity, and 

B x p y ] = ^-(VV)-±(Vn x )(n y V)~(n x Vn y ) 
2 

- -(n x Vn y V) - 2oJ c (n x T>n x T>) 



1 - lo 2 t 2 



(T>n x T>n y T>} 



2u, 



-{Vn x Vn x V) 



(B.2) 



for the Hall resistivity. 

The propagator T>((f>, <fi') in the case of white-noise dis- 
order can be expressed through the propagator 2?o(<^> <^'), 
obeying the Liouville-Boltzmann equation with only 
scattering-out term present in the collision integral, 

r i i \ 9 1 

— HjJ + iqvF cos(0 — <p q ) + (jJ c —- H — 



= 2tt<5(0 -</>'). (B.3) 

As in a zero magnetic field, the total propagator is given 
by the sum of the ladder-diagrams (thus including the 
scattering- in processes), yielding 

#1 dfo T> (<l>,<l>i)T>o(<h,<t>') 



2?r 2vr t- (D ) 



(B.4) 
(B.5) 



which we write symbolically as follows 

t - go 

Here we introduced a short-hand notation 

g (u, q) = (V ) = J g^2? o (w,q;0,^) (B.6) 

for the angle-averaged scattering-out propagator. It 
turns out that for a white-noise disorder the kernels B xx 
and B$ can be expressed in terms of go (and its deriva- 
tives with respect to q and to). The solution of HB.3(I is 
given by 

T> (u),q; </),(/)') = exp{iqR c [sin((/)' - <j> q ) - sin(<£ - <j> q )]} 

exp[wi(</> — 4>')\ 



oo 

E 

n— — oo 



—i(u) — nuj c ) + 1/r 



(B.7) 



It is worth mentioning that in the mixed-disorder model 
introduced in Sec. with both, white-noise and smooth 
disorder present, the solution of the Liouville-Boltzmann 
equation also has the form (|B.4I) . In that case, the propa- 
gator T>q satisfies the Liouville-Boltzmann equation for a 
purely smooth disorder (considered in Appendix[DJ) with 
the replacement uj — ► lu + z/r wn , where r wn is relaxation 
time due to white-noise potential. 

Using l|B.7(l and a series representation for the Bessel 
functions, we find [see, e.g. Ref. 0y| 

x , J 2 (qR c ) 

5o(^,q) = 



/! — n 

in J^qR^J-^qRc 
uj c sin7r/i 



(B. 
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where Jp(z) is the Bessel function and 



(B.9) 



In the absence of magnetic field (lo c = 0, R c = vf/^c — 
oo) the propagators T>o(u>, q; <f>, </>') and go(u), q) acquire 
a simple form 



T> (uj,q; (/>,<(/) = 



2ir5((j> - <p') 



-iw + qvF cos(0 — g ) + 1/r ' 

(B.10) 

1 1 



vV^f + (-ioj + 1/t) 2 S(u,q) 

(B.ll) 

To proceed further, we first reduce [using l|B.5|l ] the "ma- 
trix elements" appearing in IjB.ll) and l|B.2(l to the form 
containing only the propagators T>o, 



(V) = 



(DV) 



(Pq)t 

r- go' 

t 2 (VqVq) 

(r-3o) 2 ' 

T 2 {V n x ){n p V ) 



(r-9o) 2 ' 
{n x Vnp) = {n x V np) 

t- go 

(n x T>npT>) = (n x V n/3V } 

T(n x T> Q )(D np'Do) 



{Vn x V n(3 V) = 



(t - go) 2 ' 

t 2 (V n x V npV ) 

(r go) 2 
T 2 (Von x Vo)(V n p V ) 

(r-go) 3 



(B.12) 
(B.13) 
(B.14) 

(B.15) 

(B.16) 

(B.17) 



where (3 — x,y. Next, using (|B.7|) and (|B.3j) and per- 
forming the averaging over (j) q , we can express the matrix 
elements involving T>o via the propagator go- Introduc- 
ing the notation W = — iw + 1/r, we get the following 
^-averaged matrix elements, 



/t> t> \ -®9o 

{VqVo) = 

(V n x )(n x V ) = --—^.[l-Wgo] 
zq Vp 

( dgo^ 2 



(B.18) 



8vp \ dq J 



(B.19) 



/ r, v w 2 f <9 2 g 1 %> ^\ , g 
(n x V n x ) = —2 



(n x V n x V ) 



Avp \ dq 2 q dq J 2 

(B.20) 

W 9go 



(B.21) 



(n x V ) (V n x V ) 
(T> n x T> n x T> ) 



[1 - Wff ] 9g 
2QW 2 , <9<? ' 

1_ / (Pgo I dgp 

iv F \ dq 2 q dq 



f8» 

for the "longitudinal correlators" , and 



(V n x )(n y V ) 



2q^ 1 - W9 ° 



(n x V n y ) = uj c 



J F 
W dgp 
2qv 2 dq ' 



dgo 
dq ' 



(n x V nyV ) = 
(n x T> ) (T> n y T> ) = 
(V n x T> nyT> ) = 

for the "Hall correlators" . 



J F 

u c ( d 2 go , 1 dg 



Avp V dq 2 ' q 



w c ( dgo 



Av 2 F V dq ) ' 



(B.22) 

(B.23) 
(B.24) 

(B.25) 
(B.26) 
, (B.27) 
(B.28) 



J_ ( dgo iW dg 
2lo c \ dcu qvp dq J ' 



(B.29) 



Substituting Eqs. (|BTT2|l - ljBT29|) in (fBTT|) and JR2J, we 

obtain the kernels B xx and B X y averaged over <j) q , 



B x P x Hu,q) = 



r- go 
2r - g Q 
2r 2 



2 (1 - WgoY 

go 



q 2 v 2 F 



i dg 



1 dg Q 



1 



dg 



t dbJ qvpr 2 dq AvpT 3 \ dq 
(l-- 2 r 2 ) ^ 



t - go 

1 / d 2 g a 1 dg c 



22/1 1 -9o 



AvpT 2 \ dq 2 q dq 



1-OJt 



(B.30) 



T" - .90 

i / %) , iW <9g 



Auj c t 2 \ duo qvp 



2qv 2 F T 



i->v»-£ 



9g 



w c5o / 5 2 g , 1 dg 



AvpT 2 \ dq 2 q dq 



u c t + g ( dgo 



Avpr 2 r - g \ dq 



(B.31) 



In zero magnetic field, we put lo c — and substitute 
go = 1/S in l!B.30|) . After some algebra, we reduce the 
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obtained expression for the kernel B xx to the form 



(qv F ) 2 Hqv F ) 2 



2t 3 S* 3 (S- 1/r) 3 4t 2 S* 3 (S - 1/r) 2 
S-W (25- 1/t)[S- W} 2 



tS(S-1/t) 2 2T(qv F ) 2 S(S- 1/r) 2 ' 

(B.32) 

which agrees with Eq. (16b) of Ref. [ijj up to an overall 
factor l/2r related to different normalization. In the 
ballistic limit. 



Tt ^> 1, expanding i|B.32|) in 
finds the leading contribution [C(l/r)] given by the last 
two terms in (IB.32I) . 



S a + iui [Sq + ilj] 



tS$ 



So 



iuj 



rS 2 



r(qv F ) 2 S 2 ' 
1 1 
S S - i(ui + iO) 

211/2 



(B.33) 



wh ere S a = [q 2 v 2 F - (w + iQ) 2 ] 1 ' 2 . Substituting pH3|) 
in (|2.39l) and using (|C.1|I for exchange interaction, wc 
reproduce the linear-in-T correction to the resistivity in 
the ballistic regime, 



5 pZx 

Po 



T 



(B.34) 



Within the approximation of isotropic interaction used 
m Ref. El the Hartree term is determined by the triplet 
channel and is given by 



Po 



3Fg 



T 



1 + FZ E F 



(B.35) 



It is worth noting that one should exercise a certain cau- 
tion when comparing the experimental data with the 
results (IB.34|I and l|B.35|) . even in systems with short- 
range impurities. First, the higher angular harmonics 
Fn£o °f the interaction^ (neglected in the above approx- 
imation) may change the numerical coefficient in front 
of the Hartree term (see discussion in Sec. IIII El and in 
Ref. Il8j) . Second, anisotropy of the impurity scattering 
introduces an extra factor 2m>W (tt)t ^ 1 (where W(tt) 
is the effective impurity-backscattering probability) in 
both exchange and Hartree terms (see Sec. Ill Cfe and Ap- 
pendix[UJ) . The anisotropy may arise due to some amount 
of smooth disorder present in any realistic system, due to 
a finite range of scatterers, or due to the screening of orig- 
inally point-like impurities (see Sec. IIV|I. Therefore, the 
interaction parameter Fq extracted from the measured 
linear-in-T resistivity with the use of (|B.34|) , IB. 351) may 
differ considerably from that found from a measurement 
of other quantities (e.g. the resistivity correction in the 
diffusive limit or the spin susceptibility). 

To find the leading contribution to B x P y in the limit of 
vanishing magnetic field, we have to expand the propa- 
gator go up to the second order in lo c in the first term in 



curly brackets in i|B.31|) . This can be easily done by treat- 
ing the term ui c d/d(f> in l|B.3[) as a perturbation, which 
yields 



g o {B^0) = g (Lu,q;B = 0)+Lu 2 c h(Lu,q) 



1 

S 



q 2 v 2 F {S 2 - 5W 2 



8S 7 



(B.36) 



After a simple algebra, we find B X P J in the following form, 



B { x p v \u,q) 



(qv F ? 



t 2 S 3 (S - 1/r) 3 
W[S-W} 2 



(gv F ) 2 [2S - 5W] 
4t 2 S 5 {S -1/t) 2 

' 2t 2 S 4 (S-1/t) 2 ' (R37) 

which agrees with Eq. (16a) of Ref. [H In the ballistic 
limit, Tt > 1, the leading contribution [C(l/r 2 )] to B X P J 
has the form 



(j c (Sa + iu) 



Ar 2 Sl 



[<oSq - 3iS oj + 5cj 2 ]. (B.38) 



In arbitrary magnetic field, Eq. (|B.30|I can be also sig- 
nificantly simplified when the condition of the ballistic 
regime, Tt 3> 1, is assumed. Then the leading contribu- 
tion to the longitudinal MR, Ap xx — p xx (B) — p xx (0), is 
determined by the kernel 



9 2 o 



(1 - W 5 o) 2 



. dgp 

dui 



q 2 v F 



dgo 
dq 



(B.39) 



The remaining terms in (|B.30|I yield the contributions to 
the MR which are smaller at least by an additional factor 
(Tt)- 1 . Using ljB~36|) [which tells us that for uj c < T the 
magnetic-field-induced corrections to the propagator go 
are small by a factor (uj c /T) 2 ), we find that the MR for 
not very strong magnetic fields, lo c -C T, is determined 



>(p) 

>XX > 



by a quadratic in ui c correction to the kernel B x 

■ 2 ^2h(uj,q) 2h(uj,q)W(S- W) 



ABj>>( W ,g) = 



dh(u),q) _ 
dui 

oS-W 



1 



lv 2 F S* 



4r5 6 



S 2 



q 2 v 2 F 



dS_ 

dq 

5W 2 



S 



S 



5W{S + W)(3S 2 -7W 2 ) 



-0J r 



,Sq 



2S 3 
iui 



S + W 



8tS$ 



45^ + 13iuS% 



25uj 2 So - Sbiuj^So + 35w 4 



(B.40) 



independently of the relation between ui c and r 1 . Sim- 
ilarly, using i|B.36|) . one can find the correction to 
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Eq. 



38(1 in a finite magnetic field oj c <C T, 

, /LA. I.' - ' -'- 1 
(/,, ^ — _/,, J _ 

y 



A5W( Wl9 ) = -u f^l+^ pSZ + 7c 2 ]. (B.41) 

Again, this correction is independent of the relation be- 

-l 



tween u> c and r . The results (IB.40|I and (|B.4ip are used 
for calculation of the interaction-induced corrections to 
p xx and p xy for the white-noise disorder and oj c <C T in 
Sec.lVBl 



APPENDIX C: LINEAR-IN-T TERM IN THE 
BALLISTIC LIMIT AT B = 

In this appendix we calculate the leading ballistic cor- 
rection to the conductivity at B = for a generic scat- 
tering cross-section W(<f>— <f>') in the case of the Coulomb 
interaction. As explained in Sec. Ill CR . this term (pro- 
portional to Tt) is obtained by substituting the bal- 
listic asymptotics (I2.51|) of B xx in the general formula 
((2.35|) . Likewise, the interaction propagator U(u>,q) en- 
tering ((2.35|) has to be replaced by 



U(u,q) 



2v l + iuj(Vi) 



(C.l) 



with the free propagator T>[ given by Eq. ((2.50(1 . Per- 
forming the angular integration (...), we get 



~ 2ir 2 
QdQ 



Tt Im 



d ( ft 
dft-— ftcoth— 

on \ 2 



^-(ft + zO) 2 ] 1 ^ 



[Q 2 - (ft + iO) 2 ] 1 / 2 + i(ft + iO) 
d(j) (-iQ)W((f>)(l-cos(t>) 
2^ [Q 2 cos 2 1 - (n + M) 2 }[Q 2 - (ft + iO) 2 } 1 / 2 

-in 



[Q 2 - (ft + iO) 2 ] 3 / 2 



(C.2) 



where we introduced the dimensionlcss variables ft — 
u/T, Q = qv F /T, and W((f>) = 2ttvtW {<f>) . It is con- 
venient to split the interaction propagator as follows 



2vU(Sl,Q) 



Sq 



So + i(Q + iO) 



1 



s a + i(n + io) 

(C.3) 

where S = [Q 2 - (ft + z0) 2 ] 1/2 . The first term corre- 
sponds to a statically screened interaction and is equiva- 
lent to a point-like interaction with Vb = 1/2;;, the second 
term results from the dynamical weakening of screening. 

As discussed in Sec. Ill CB . the contribution Sai 1 ] of the 
first (constant) term is proportional to the backscatter- 
ing probability W(tt), see Eq. ((2.52(1 . Let us show that 
this follows also from Eq. ((C.2(l . Performing the variable 
change Q — > So, we find 



- ~^ Tt Im 



dft— ( ftcoth- ) $(ft) 







-Q. 



ImSo 



ReSo 



C 



i C 
t 



FIG. 14: The contours C and C' of integration over So in 
Eq. SUll 



Hft)(s 2 + ft 2 ) 



$(ft) = / dS 



d. ( p 



o2 
°0 



2sin 4 ((/)/2)I4>(0) 



2tt S 2 cos 2 f - (ft 



iO) 2 sin 2 | ' 



(C.4) 



The contour C of integration over So in Eq. I(C.4() is 
shown in Fig. 1141 Interchanging the order of integration 
over (j) and So, we see that for any cf> n (i.e. cos | ^= 0) 
the So-integral converges. Furthermore, transforming the 
integration contour C — > C as shown in Fig. it is 
straightforward to reduce $(ft) to an explicitly real form. 
Therefore, only the singular point <f> — tt (where the result 
of So-integration diverges as 1/| cos % |, implying that the 
imaginary part of $(ft) is determined by a delta-function 
in /^-integral) contributes to l(C.4() . so that Sai 1 ] cx W(ir). 
To find the corresponding coefficient, one can consider 
the isotropic scattering W(4>) = const and to integrate 
over cj) first, yielding 



W(tt)Tt, 



(C.5) 



in agreement with ((2.52(1 . Note that the integral over 
ft is formally divergent at the upper limit. It should be 
cut off at ft ~ EpjT yielding a temperature independent 
contribution ~ e 2 W(rr)EpT which renormalizes the value 
of the Drude conductivity. 

(2) 

We now turn to the contribution oo xx of the second 
(dynamical) term in the interaction propagator ((C.3(l . 
which differs from Eq. I(C.4I) by an extra factor — zft/[So + 
i(n + iO)]. Rotating at <f> ^ tt the integration contour as 
before, we reduce the So-integral to the form (So — > —iY) 



dY 



Y + ft 



Y 2 (Y 2 cos 2 | + ft 2 sin 2 |) : 



(C.6) 



which is again real and thus yields no contribution to 

(2) 

5<J XX ■ Though the point <f> — tt is singular in this case as 



.32 



well, the singularity is only logarithmic (~ In | cos ^|), so 
that no contribution proportional to W(n) arises. This 
can be easily checked by assuming W(4>) = const and 
performing the 4>- integration first. Therefore 

s4 2 J = o, 

and the linear-in-T term is given by Eq. JHSJ. 

In the above consideration we have expanded the bal- 
listic propagator T> up to terms with one scattering event. 
In the case of small- angle scattering this is justified pro- 
vided Tt s 3> 1, while in the intermediate temperature 
range r _1 <C T -C t~ 1 processes with many scatter- 
ing events dominate (though the particle motion is typ- 
ically close to the straight line). The term Sa xx which 
is governed by anomalous processes of returns in a time 
t<T _1 <Cris exponentially small in this case, see 

Sec. lII Ck . As to the Saxx contribution to the linear-in-T 
term, it remains zero in this case as well. To demonstrate 
this, we use Eq. (|2.36l) . In the first and the third terms we 
can replace T> by the free propagator (|2.5U|) . the fourth 
term gives no Tt contribution, while in the second term 
we should take into account the angular diffusion (|2.9|l 
around the straight trajectory, 



1 



-<£>) - (n x Vn x ) 



2r[q 2 v 2 F - (u + i0) 2 ] 3 / 2 



(C.7) 



Combining the contributions to B xx of all the three 
terms, we get 



B xx (uj,q) 



—ioj 



[q 2 v 2 F - (w + i0) 2 } 3 / 2 \2 2 



1 1 



l)=0, 
(C.8) 



so that the coefficient of the Tr-term indeed vanishes. 



APPENDIX D: SOLUTION OF 
LIOUVILLE-BOLTZMANN EQUATION FOR A 
SMOOTH DISORDER 

In this appendix, we will solve the classical equation 
for a propagator of a particle moving in a smooth random 
potential in a magnetic field, 



d 



1 d 2 



iqVF COS ( 



which is just a standard expansion in eigenfunctions of 
the Liouvillc-Boltzmann operator. We now treat the 
term 5C — ivpq cos <j> as a perturbation. The first-order 
correction to the eigenvalues A^ = — iuj + inuj c + n 2 /t 
vanishes, while the second order correction is 



Ai 2 



2, ,2 



ft) 



1 



2 1 - (iuj c t + 2n) 2 



(D.3) 



and can be neglected along with — iui in all terms except 
for n — in the diffusive limit. The first order correction 
to the right eigenfunction for n — reads 



iqvpT 
1 + uj 2 t 2 



(D.4) 



while the left eigenfunction differs from l|D.4(l by a re- 
placement uj c — > —lu c . Thus, in the diffusive limit the 
propagator has the form 



2?(w,q;<M') 



1 



Dq 2 — iuj 



1 - 



i qvF t (cos 4> + w c t sin <j> ) 



1 + UJ 2 T 2 

^ i qvF t (cos (j)' — lu c t sin <f>') 

1 + UJ 2 T 2 



E 



e in{4>-4>') 

inuj c + n 2 It 



(D.5) 



In a strong magnetic field (uj c t 1) one can go be- 
yond the diffusion approximation. In this case one can 
represent the propagator in the form 

V{oj, q: 0, 4>') — d(u>, q; <j>, (f)') exp[—iqR c (sin(f) — sin <£')], 

(D.6) 

and solve the equation for d(4>, (/>'), 



.qv F . 
-iuj — i sin < 



1 V F 

uj r + 2i cos ( 



d 



1 / qv F 



i d 2 



cos 



Td<j) 2 



(D.7) 



perturbatively in q. At q — we have the same solution 
(|D.2() as in the diffusive limit. The first order correc- 
tion to the eigenvalues is now produced by the (7 2 -term 



= 2wS((f) - <f>'). (D.l) in (ITXTI) . \ { n ] = Dq 2 , with D = R 2 c /2t the diffusion 



Here the polar angle of the velocity is counted from the 
angle of q, <f) — <fiq — * <fi- 

We first consider the diffusive limit, Dq 2 , lo <C 1/t, 
and solve this equation perturbatively in q for arbitrary 
magnetic field. Putting q = 0, we obtain the solution in 
the form 



constant in a strong magnetic field. The second order 
corrections turn out to be small compared to 
for (qR c ) 2 <§C u c t. As in the diffusive limit, for calcula- 
tion of B xx the corrections to the eigenfunctions ^ ra with 
n 7^ can be neglected. The first-order correction to 
is found to be (we drop the term cx sin 20, since it does 
not contribute to B xx in the leading order) 
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in(4>-4>') 



-iuj + inuj c + n 2 jr 



(D.2) 



tqvFT cos ( 

(^cT) 2 



(D. 



33 



leading to Eq. (|3.1|l . 

To calculate the kernel B xy , we need a more accurate 
form of the propagator. Therefore, we should analyze 
the corrections to the eigenvalues and eigenfunctions of 
the Liouvillc-Boltzmann operator to the next order in 
(qR c ) 2 /uj c t. To do this, it is convenient to perform the 
transformation 



£>(w,q; </>,</>') = d(cj, q; 6, </)') 

. qRc 



x exp 



'1 + /3 2 
+ /3(cos 6 — cos 0')]} 



[/3 2 (sin0-sin0') 



(D.9) 



and introduce the dimcnsionless variables (3 = lo c t, Q — 
qR c (3/(l + (3 2 ) 1/2 , n = 2ujt. The equation for 
d(u, q; 6, 6') takes then the form 



-i— + Q 2 cos 2 . 



+2iQcos 



t d 



86 
d 2 



<w,q;0,^') 



2-kt8{6- 6'), 
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where we performed a rotation 6 = 6+6p, 6p = arccot/3. 
Treating for Q 2 <C max[l, (3} (i.e. [qR c ) 2 -C uj c t in a 
strong magnetic field, (3^-1) the term 



Q 2 ~ ~ 8 

SL = — cos 26 + 2iQ cos 6 — = 
2 86 

as a perturbation to the operator 
.fl Q 2 



(D.ll) 



L 



2 2 

we find the unperturbed solution 



(D.12) 
(D.13) 

and the first-order correction to the eigenvalues An = 0. 
Calculating the n — eigenfunctions and eigenvalue up 
to the second order in the perturbation (|D.11J| we finally 
obtain the singular part of the propagator for f3 1 with 
required accuracy, 



-itt + Q 2 + 2in/3 + 2n 2 



V s (ui, q; 6,6') = 2rcxp[— iQ(sin6 — sin 6')] 
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where Q = qR c and the functions \R l (0; Q) are given 
by 



Q 2 

i<5 cos ± — sin 20 
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— ±iQ{l + — )sin( 



5Q 2 7iQ 3 
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4 r 24 ^ 
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64 



(D.15) 



As to the regular part of the propagator, for n ^ it is 
sufficient to calculate the eigenfunctions to the first order 
in the perturbation, which yields 



X> rog (w,q;0,0') = 2rexp[-iQ(sin0-sin0')] 
-in + Q 2 + 2inf3 + 2n 2 



(D.16) 



where 



R,L\ 



?;n) = 1 - — cos 

T — sm20± — — sin0. 
4/5 p 

(D.17) 

The results (|D.14f) - (|D.17f) allow us to calculate the ker- 
nel B xy (u> , q) in the first non- vanishing order in (3~ l , see 
Sec. Pirn! 



APPENDIX E: PROPAGATOR FOR 
ANISOTROPIC SYSTEMS 



In this Appendix, we assume that the collision integral 
C induces a transport anisotropy, i.e. that the scattering 
cross-section W(6, 6') is not a function of 6 — 6'. The 
propagator T>(u>, q; 6, 6') satisfies the equation 



-ico + iqvF cos( 
= 2w6(6-6'), 
where 

[£*](0) = i/ 



, \ 9 a, 

6 q ) + WcTTI + ^ 



2?(w,q; 0, 0') 
(E.l) 



2tt 



[*(6)-*(<p')]W(6,6'). (E.2) 



We first consider the diffusive limit and concentrate on 
the leading contribution T> s governed by the diffusion 
mode. 

This requires finding a lowest eigenvalue A of the op- 
erator in the l.h.s. of (|E.1|I and the corresponding left and 
right eigenfunctions. Treating the term iqvF cos(0 — 6') 
perturbatively as in Appendix El we find 

^r,l{ 6 ) = 1 - W>f ( ±UJ c^- + 6) cos(0 - 6 q ) (E.3) 



and 



A = D a f}q a qf3 - iuj, 



(E.4) 
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with the diffusion tensor 



(E.5) 



According to l|3.3|l - i|3.5[l . we on ly need to calculate aver- 
ages of the type (D) and (n a T>), so that it is sufficient 
to keep the zero and first harmonics in 4> in Eq. I jfi.lOjl . 
Using 



We thus get the result l|6.5|l for the singular contribution 
V s , with ty RtL given by l|E.3|) . 

In a strong magnetic field (lo c t ^> 1), we can go be- 
yond the diffusive limit. Proceeding as for an isotropic 
system, we perform the transformation l|D.6|) . Treat- 
ing the q-dependent terms in the obtained equation for 
d{oj, q; </>, (j)') as a perturbation and keeping the singular 
contribution only, we come to the result (|6.1U|) . where 
x{4>) can be represented symbolically as 



(n a Cnp) 



(E.7) 



x(<P) 



iqvF ( d 



Csin 



(E.6) 



we then reduce (|E.6(I to the form l|6.11(l . 
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